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STUDY  OF  BRIGHTNESS  AND  CURRENT  LIMITATIONS  IN  INTENSE 

CHARGED  PARTICLE  BEAMS 


M.  Reiser  and  S.  Guhairay 

University  of  Maryland,  College  Park,  MD  20742  (301)  405-4960 
INDIVIDUAL  RESEARCH  ABSTRACT 

Over  the  past  several  years  of  ONR  support  for  our  research  program  we  have  mainly 
studied  the  various  schemes  for  intense,  high-brightness  H~  beam  transport  and  focusing  in 
the  context  of  its  application  in  space  defense.  Detailed  theoretical  studies  revealed  that  the 
conventional  gas  focusing  system  is  not  suitable  as  a  low-energy  beam  transport  (LEBT) 
system  and  also  that  there  are  too  many  unknown  parameters  to  model  accurately  the  be¬ 
havior  of  partially  charge-neutralized  particle  beams.  We  concluded  that  the  electrostatic 
quadrupole  lens  system  will  be  a  good  choice.  We  have  developed  a  large  number  of  sim¬ 
ulation  codes  and  also  accessed  into  the  existing  codes  in  the  accelerator  community  (e.g., 
PARMILA,  SNOW-2D,  PARMTEQ,  etc.)  to  strengthen  our  analysis. 

During  the  1992-93  contract  period  we  focused  our  attention  to  the  experimental  activities 
on  H"  beam  characterization  and  on  the  installation  of  a  LEBT  system  for  beam  transport 
experiments.  We  have  simultaneously  improved  our  code  by  incorporating  many  practical 
features  that  we  encountered  during  the  analysis  of  experimental  data. 

We  have  studied  H~  beams  from  two  types  of  ion  sources:  a  volume  ionization  type 
and  a  magnetron  type  source.  One  of  the  major  problems  in  this  work  is  to  transform  a 
highly  diverging  beam  from  the  source  into  a  highly  converging  one  so  that  the  output  beam 
from  the  LEBT  can  be  matched  into  the  acceptance  ellipse  of  an  RFQ.  Furthermore,  the 
emittance  budget  is  quite  restricted.  After  a  detailed  beam  dynamics  studies,  we  have  been 
finally  successful  to  design  an  optimized  LEBT  system  that  can  deliver  a  matched  beam  to 
the  RFQ;  this  is  the  first  time  that  such  good  beam  matching  has  been  simulated  in  the 
context  of  the  particular  experiment. 

We  have  set  up  an  experimental  facility  in-house  to  perform  various  tests  of  the  LEBT 
system.  We  have  made  several  important  modifications  of  our  ESQ  LEBT  system  which  was 
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developed  Icist  year.  The  electrical  layout  of  the  lens  assembly  has  been  changed  and  this 
required  to  modify  the  vacuum  vessel.  The  system  is  expected  to  have  a  long-term  reliability 
in  terms  of  voltage  holding.  The  ESQ  LEBT  is  supported  in  the  vacuum  vessel  using  a 
vacuum  manipulator  developed  in-house.  This  allows  to  align  the  system  in-situ.  Various 
tests  of  the  ESQ  LEBT  are  being  conducted. 

We  have  received  a  magnetron  ion  source  from  the  Superconducting  Super  Collider  Lab¬ 
oratory  (SSCL).  This  allows  us  to  set  up  an  in-house  test  stand  facility  for  H~  beams. 

In  theoretical  work  supported  by  this  contract,  Christopher  Allen,  a  graduate  research 
assistant  working  towards  his  Ph.D.  has  developed  a  new  moment-method  Laplace  solver 
for  low-energy  beam  transport  codes.  This  technique  is  expected  to  increase  both  speed  and 
accuracy  of  electrostatic  field  calculations. 
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STUDY  OF  BRIGHTNESS  AND  CURRENT  LIMITATIONS  IN  INTENSE 

CHARGED  PARTICLE  BEAMS 


M.  Reiser  and  S.  Guharay,  University  of  Maryland 


RESEARCH  OBJECTIVES: 

Intense,  high- brightness  neutral  beams  have  a  strong  potential  to  strengthen  the  space 
defense  program.  In  recent  years  a  strong  emphasis  has  been  laid  on  pushing  the  state-of- 
the-art  of  negative  ion  beams,  especially  the  H“  beams.  This  has  occurred  due  to  several 
special  features  of  the  negative  ions,  e.g.,  non-Liouvillean  stacking  to  generate  intense  beams, 
efficient  charge-exchange  process  to  generate  high-energy  neutral  beams,  and  furthermore, 
future  application  of  H~  beams  in  ion-beam  microlithography.  Such  a  demand  from  multiple 
disciplines  has  set  a  strong  challenge  to  the  community  involved  in  research  on  ion  sources 
and  ion-beam  dynamics,  and  a  rapid  progress  of  the  field  is  visible.  Our  research  program  on 
nigh-brightness  H”  beams  involves  both  experimental  and  theoretical  work  on  the  generation 
and  transport  of  the  beams;  the  major  emphasis  here  is  on  the  study  of  the  beam  emittance. 

After  a  detailed  investigation  of  the  various  methods  of  transporting  a  high-brightness 
H~  beam  we  have  determined  that  an  electrostatic  quadrupole  lens  system  is  a  good  choice. 
We  have  set  up  a  strong  infrastructure  on  simulation  of  beam  dynamics.  This  work  is  being 
continuously  examined  in  the  light  of  experimental  results.  Consequently,  many  practical 
features  are  included  in  our  simulation  codes  and  the  simulation  predictions  are  expected  to 
be  very  reliable.  Experiments  with  an  ESQ  LEBT  system,  which  was  developed  in-house,  are 
in  progress.  We  have  also  initiated  an  experimental  program  to  set  up  a  test  stand  facility 
at  Maryland. 
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STUDY  OF  BRIGHTNESS  AND  CURRENT  LIMITATIONS  IN  INTENSE 

CHARGED  PARTICLE  BEAMS 


M.  Reiser  and  S.  Guharay,  University  of  Mciryland 

STATUS  OF  RESEARCH  IN  PROGRESS 

Our  research  during  the  current  contract  period  has  been  mciinly  motivated  towards  the 
experimental  work  on  transport  of  high-brightness  H“  beams  and  focusing.  The  experi¬ 
mental  studies  are  being  carried  out  at  the  test  stand  facility  of  the  Superconducting  Super 
Collider  Laboratory  (SSCL)  using  ion  beams  from  a  volume  ionization  source:  a  new  in-house 
program  is  also  established  with  some  support  for  equipment  from  SSCL  to  set  up  a  test 
stand  using  a  magnetron-type  ion  source.  As  the  characteristic  parameters  of  H“  beams  from 
these  two  different  ion  sources  are  quite  different,  these  studies  allow  us  to  address  a  number 
of  key  physics  issues  on  beam  dynamics  as  well  as  to  evaluate  some  practical  problems  in 
terms  of  long-term  stable  operation,  ease  of  handling,  etc.  This  work  is  described  below. 

High-Brightness  H“  Beam  Transport  and  Focusing 

As  was  mentioned  earlier  in  our  previous  reports  to  ONR.  the  study  of  high- brightness 
H“  beams  is  motivated  towards  the  development  of  an  efficient  low-energy  beam  transport 
channel.  After  a  detailed  investigation  of  the  characteristics  of  various  types  of  transport 
systems,  we  decided  to  focus  our  attention  on  an  electrostatic  system,  mainly  the  electro¬ 
static  quadrupole  configuration,  in  the  context  of  transporting  short-pulse  beams  which  are 
particularly  relevant  to  our  beam  transport  experiments  with  volume  and  magnetron  sources 
(typically,  pulse  length  beam  neutralization  time).  Our  approach  is  to  initially  develop 
a  good  understanding  of  the  beam  dynamics  through  simulation  studies,  and  to  build  up 
a  close  link  between  experimental  results  and  simulation  predictions.  This  procedure  helps 
to  continuously  improve  upon  any  weak  points  of  the  simulation  schemes,  and  it,  therefore, 
enhances  the  confidence  level  of  the  simulation  predictions. 
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The  goal  of  this  work  is  to  deliver  a  matched  beam  to  the  RFQ  in  the  linac  section 
without  any  significant  emittance  dilution.  In  the  case  of  the  SSC  RFQ  the  matching  condi¬ 
tion  is  given  by  the  acceptance  ellipse  parameters:  o  =  1.26,  ^  =  1.86  cm/rad,  =  0.2t 
mm-mrad.  This  demands  a  strong  convergence  of  the  beam  from  the  LEBT  section.  We 
find  that  a  short  einzel  lens  module  between  the  ESQ  LEBT  and  the  RFQ  will  allow  us  to 
handle  the  emittance  budget  efficiently,  especially  when  a  large  drift  space  (~  few  cm)  exists 
to  accommodate  the  mechanical  structures  as  the  front  wall  of  the  RFQ;  however,  this  space 
should  be  kept  to  a  minimum  to  deliver  a  good  beam  to  the  RFQ.  It  is  estimated  using  the 
SN0W-2D  code  that  the  einzel  lens  geometry,  shown  schematically  in  Fig.  la,  can  satisfy 
the  nominal  matching  condition  at  a  distance  of  about  3.0  cm  from  the  front  end  of  the 
RFQ  wall  without  any  noticeable  emittance  growth  (Fig.  lb);  the  third  electrode  at  ground 
potential  simulates  the  front  wall  of  the  RFQ  here.  This  result  corresponds  to  the  follow¬ 
ing  beam  parameters  at  the  input  of  the  einzel  lens:  beam  current  =  30  mA,  e„  =  0.23;r 
mm-mrad,  X  ~  7  mm  and  X'  ~  —50  mrad.  Our  ESQ  LEBT  is  chosen  to  deliver  an  output 
beam  so  that  it  can  satisfy  the  aforementioned  requirements  for  matching  with  the  einzel  lens. 

A.  H“  beams  from  a  volume  source 

The  emittance  meeisurements  of  a  30  mA,  35  kV  H“  beam  from  the  SSCL  volume  source 
are  shown  as  contour  plots  in  Fig.2.  These  data  are  taken  at  a  distance  of  about  10  cm  from 
the  extraction  electrode  (aperture  radius  =  4  mm).  The  beam  parameters  at  this  location 
are:  beam  size  =  2.38  cm,  full  divergence  A6  =  260  mrad,  7re„  =  0.15377r  mm-mrad.  The 
drift  space  of  10  cm  long  is  due  to  the  insertion  of  an  electron  suppression  system  which 
separates  out  the  electron  component  (initial  electron- to- ion  current  ratio  is  about  40.)  from 
the  extracted  current.  Such  a  long  drift  space  causes  a  significant  blow-up  of  the  H"  beam 
envelope.  It  is  currently  planned  to  replace  the  electron  suppression  system  by  a  shorter  one 
of  about  5  cm  long. 

In  order  to  estimate  the  beam  parameters  at  the  extraction  electrode  the  aforementioned 
beam  data  at  z  =  10  cm  are  given  as  input  to  the  envelope  simulation  code  which  solves  the 
K-V  equations.  This  code,  developed  here,  has  a  number  of  important  features  that  allow 
to  include  several  practical  considerations,  e.g.,  the  space-charge  effects  due  to  the  accum- 
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accumulation  of  extraneous  charge  elements  (electrons  and  positive  ions)  in  the  neighborhood 
of  the  extraction  region,  plausible  profiles  of  the  extraneous  charge  elements,  etc.  Figure  3(a) 
shows  the  eissumed  space-charge  correction  factor,  /,  due  to  the  electrons.  Note  that  /  is 
negative  here,  and  the  beam  perveance  is  to  be  multiplied  by  a  factor  of  ( 1  —  /)  to  include  the 
effect  of  the  electrons  in  the  space-charge  force  term  in  K-V  equations.  The  estimated  beam 
envelope  upstream  towards  the  extraction  electrode  (here,  at  2  =  0)  is  shown  in  Fig.  3(b). 
This  suggests  that  the  beam  forms  a  waist  close  to  the  extraction  electrode  and  it  fills  up 
almost  the  full  aperture. 

Using  the  above  information  of  the  beam  parameters  at  the  extraction  electrode,  the 
characteristics  of  the  ESQ  LEBT  system  and  the  behavior  of  the  beam  through  the  LEBT 
are  determined.  Figure  4  shows  the  envelope  of  the  H~  beam  through  the  ESQ  LEBT.  It 
may  be  noted  that  the  aperture  of  the  ESQ  lenses  are  about  a  factor  of  two  larger  than  what 
we  mentioned  in  our  last  year’s  report;  this  change  has  been  necessary  to  accommodate 
the  large  beam  at  the  end  of  the  5  cm-long  electron  suppressor  in  this  experiment.  The 
distribution  of  the  beam  particles  through  the  ESQ  LEBT  is  estimated  using  the  modified 
PARMILA  code.  Figure  5  shows  the  particle  distribution  in  phase  space.  The  output  beam 
parameters  are:  X  =  7.2  mm,  Y  =  7.3  mm,  X'  =  —51  mrad,  Y'  =  —51  mrad,  injjtni  ~  1.5. 
The  emittance  growth  is  mainly  due  to  chromatic  aberrations.  Note  that  the  output  beam 
parameters  match  closely  to  the  desired  characteristics  of  the  input  beam  for  the  einzel  lens. 

Details  of  this  work  are  given  in  reference  4.  Note  that  the  reference  numbers  referred 
to  in  this  section  correspond  to  the  numbers  in  the  list  of  presentations,  publications  and 
invited  talks  (p.  23)  of  this  report. 

B.  H~  beams  from  a  magnetron  source 

The  emittance  measurement  of  the  H~  beam  from  the  SSCL  magnetron  source  is  shown  in 
Fig.  6.  These  data  were  taken  at  11.75  cm  downstream  from  the  extraction  slit  of  the  source; 
the  beam  size  is  about  3.5  cm,  full  beam  divergence  is  about  300  mrad,  and  7re„  =  0.12x 
mm-mrad.  In  the  case  of  a  magnetron  source  the  ratio  of  e/H~  in  the  extracted  current  is 
about  1-2;  hence  the  space-charge  effect  of  the  electrons  may  not  be  important.  Following  the 
method  as  described  in  the  previous  case,  the  beam  parameters  at  the  extractor  is  estimated 
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using  the  aforementioned  emittance  data,  and  it  is  found  that  X  ~  1.1  mm  and  X'  ~  70 
mrad. 

The  beam  envelope  through  the  ESQ  LEBT  is  shown  in  Fig.  7.  It  is  noted  from  the 
particle  simulation  results  that  the  emittance  growth  for  the  full  beam  current  is  quite  large 
here  -  a  factor  of  about  3.  Scraping  off  about  15%  of  the  beam  particles  the  emittance 
growth  drops  to  a  factor  of  1.5;  Fig.8  shows  the  phase-space  distribution  of  the  particles. 
The  beam  parameters  are:  X  =  4.7  mm,  Y  =  4.4  mm,  X'  =  —49  mrad,  Y'  —  —48  mrad. 
These  results  also  match  reasonably  well  with  the  requirements  of  the  einzel  lens  input  as 
discussed  earlier. 

Details  of  this  work  are  given  in  reference  1. 

C.  ESQ  LEBT  system  and  tests 

It  is  planned  to  initially  use  a  prototype  ESQ  LEBT  system  which  has  already  been  de¬ 
veloped  in-house  at  Maryland.  The  compatibility  of  the  LEBT  system  with  the  SSC  volume 
source  environment  has  been  checked.  Figure  9  shows  an  initial  experimental  set-up  of  the 
LEBT  on  the  SSCL  test  stand.  Figure  10  shows  the  test  facility  at  Maryland  to  perform 
voltage  hold-off  tests.  The  experiment  is  now  in  progress. 

D.  Magnetron  ion  source  facility  at  the  University  of  Maryland 

We  have  received  a  magnetron  ion  source  with  numerous  supporting  electronics  from 
SSCL.  Figure  11  shows  a  few  major  components  of  the  experimental  facility.  In  order  to 
bring  the  test  stand  to  operation  and  run  the  ion  source  meaningfully,  it  will  be  necessary 
to  develop  the  computer  control  system,  and  also  some  diagnostic  equipment,  e.g.,  Faraday 
cup,  toroids,  emittance  scanner,  etc.,  should  be  procured. 

E.  Theoretical  studies 
1.  Laplace/Poisson  Solver 

An  efficient  Laplace/ Poisson  Solver  has  been  developed.  It  is  fully  3D  and  capable  of 
running  practical  problems  on  an  IBM  PC.  The  technique  relies  on  a  combination  of  the 
method  of  moments  and  fast  iterative  techniques  for  solving  linear  systems.  Specifically,  we 
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proceed  by  reformulating  Laplace’s  equation  into  an  integral  equation  over  the  boundary 
surface,  reducing  the  dimensionality  of  the  original  system.  The  new  problem  is  approxi¬ 
mated  by  the  method  of  moments  to  yield  the  matrix-vector  equation  =  y.  VVe  then  use 
conjugate  gradient  algorithms  to  solve  this  equation;  this  is  an  iterative  method  which  seems 
to  provide  fastest  convergence. 

The  above  technique  haa  been  implemented  on  an  i486  PC  operating  at  33  MIL  and 
running  Windows  3.1  operating  system.  A  3D  potential  problem  is  the  modeling  of  an 
electrostatic  quadrupole  lens.  Figure  12  shows  the  computer  model  of  an  ESQ  lens.  It  is 
formed  from  4  cigar-shaped  electrodes,  the  beam  would  enter  from  the  left.  Each  electrode 
is  59  mm  long  and  has  a  radius  of  12  mm.  The  aperture  of  the  entire  lens  is  10.5  mm.  Two 
grounding  shunts  are  located  at  r  =  ±  31  mm  (they  are  not  shown  in  Figure  12  to  avoid 
clutter)  which  provide  isolation  from  adjacent  lenses. 

The  single  particle  focusing  effect  (the  kappa  function  ^(2))  from  such  lenses  can  be 
determined  from  the  on  axis  derivatives  dExfdx  and  dEy/dy.  Figure  13  shows  the  computed 
data  for  the  y-plane  for  the  case  in  which  the  i-plane  electrodes  are  driven  to  1  V  and  the 
y-plane  electrodes  are  held  at  —1  V.  The  grounding  shunts  at  either  end  of  the  lens  cause 
the  rapid  decay  in  dEyjdy. 

We  can  use  the  Poisson  e.xtension  to  simulate  a  uniform  charge  density  ellipsoid  in  a 
conducting  cylinder.  This  situation  is  useful  in  modeling  cold  bunched  beams  propagating 
through  a  beam  pipe.  There  exists  an  analytic  solution  for  the  potential  of  such  an  ellipsoid 
in  free  space.  Thus,  it  is  only  necessary  to  model  the  pipe  (surface  charge)  numerically.  The 
above  simulations  can  be  used  to  determine  the  so  called  “g-factor”  for  bunched  beams  in 
cylindrical  pipes. 

Details  of  this  work  are  given  in  reference  3. 
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Figure  1(b).  SNOW-2D  results  of  the  particle  distribution:  input 
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Figure  2.  Contour  plots  of  the  beam  from  the  SSC  volume  source  at  a  distance  of  10  cm 
downstream  from  the  extractor. 
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Figure  3.  Estimation  of  beam  characteristics  at  the  extractor. 
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Figure  4.  K-V  envelope  solution  for  an  H  beam  from  the  volume  source. 
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Figure  5.  (a)  Particle  distribution  at  the  input  of  the  ESQ  LEBT;  (b)  particle  distribution 
at  the  output  of  the  ESQ  LEBT. 
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Figure  6.  Contour  plots  of  the  beam  from  the  SSC  magnetron  source. 
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Figure  8.  Particle  distribution  from  modified  PARMILA. 
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ABSTRACT 

-A  compact  6-lens  electrostatic  quadrupole  lens  system  in  conjunction  with  a 
short,  single  einzei  lens  section  has  been  developed  with  the  aim  of  transporting  a 
30  mA,  35  kV  H~  beam  (normalized  beam  brightness  of  about  10"  A/(m-rad)*) 
over  a  length  of  about  30  cm  and  focusing  it  into  an  RFQ.  The  effect  of  a  neutrad 
background  gcis  on  the  meaisurements  of  the  beam  from  the  ion  source  is  studied 
in  order  to  evaluate  reliably  the  input  beam  parameters  for  the  lens  design.  The 
beam  dynamics  calculations  have  been  made  using  simulation  codes  for  typical  H“ 
beam  parameters  from  two  different  ion  sources:  a  magnetron  type  and  a  volume 
source  type.  The  simulation  results  show  a  relatively  modest  emittance  growth, 
that  is  within  about  50-60%  even  in  the  case  of  an  input  beam  with  very  small 
beam  radius  (1.1  mm)  and  very  large  divergence  (72  mrad). 

INTRODUCTION 

Over  the  past  few  years  our  research  hcis  been  primarily  focused  on  simulation 
studies  of  charged-particle  beam  dynamics  with  the  aim  of  developing  an  effi¬ 
cient  low-energy  beam  transport  (LEBT)  system  for  intense,  high-brightness  H~ 
beams.  While  dealing  with  the  problem  of  design  and  development  of  an  LEBT. 
we  encountered  several  problems  of  practical  interest.  One  of  the  puzzling  issues 
is  related  to  the  assumption  of  input  beam  parameters.  Often,  the  beam  param¬ 
eters  are  measured  at  a  certain  distance,  typically  ~10  cm.  downstream  from  the 
extraction  aperture  of  the  .  source.  This  leads  to  some  uncertainty  in  predicting 
beam  parameters  at  the  extraction  point  as  this  information  is  required  to  define 
the  input  beam  parameters  in  designing  a  LEBT.  Second,  tne  beam  parameters 
at  the  e.xtraction  aperture  depend  very  much  on  the  mechanism  of  ion  source 
operation  and  the  extraction  optics.  It  is  thus  of  practical  interest  to  adapt  an 
LEBT  scheme  which  may  function  reasonably  well  for  a  certain  range  of  the  input 
beam  parameters.  Such  a  system  is  less  susceptible  to  failure  in  th»,  event  the 
source  parameters  drift  to  a  certain  extent;  this  situation  is  encountered  quite 
commonly  in  experiments  with  ion  sources.  Third,  in  the  injectors  of  high-energy 
accelerators  a  good  buffer  space.  ~20  -  .30  cm.  between  the  ion  source  and  the  first 
stage  of  acceleration,  commonly  an  RFQ  (radio-frequency  quadrupole),  is  desired 
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as  it  offers  some  practical  advantages  in  terms  of  clean  operation  of  the  RFQ  and 
flexibility  in  achieving  satisfactory  transformation  of  beam  parameters. 

The  present  article  addresses  the  aforementioned  problems  when  widely  dif¬ 
ferent  beam  parameters  are  used  as  input  conditions.  Special  considerations  are 
given  to  meet  the  matching  condition  dictated  by  the  acceptance  ellipse  parame¬ 
ters  of  the  RFQ  in  the  Superconducting  Super  Collider  (SSC).‘ 

INPUT  BEAM  PARAMETERS  AND  THE  ROLE  OF  NEUTRALIZING 

BACKGROUND 

Beam  characterization  experiments  were  done  on  the  SSC  test  stand  with  H“ 
beams  from  both  magnetron  and  volume  ion  sources.  In  this  article,  we  will 
mainly  deal  with  the  magnetron  beam  parameters;  some  analysis  is  done  using 
typical  beam  parameters  of  the  Brookhaven  National  Laboratory  (BNL)  volume 
source.  Sadaatmand  et  al.  reported  some  recent  results  on  the  characteristics 
of  H~  beams  from  the  SSCL  volume  source*;  we  plan  to  use  these  results  in  the 
future  to  study  beam  dynamics  in  the  context  of  our  LEBT  system. 

Figure  i  shows  the  pha^e-space  distribution  of  the  beam  particles,  when  the 
emittance  scanner  slit  was  located  at  11.75  cm  downstream  from  the  extraction 
slit  of  the  source.  Using  the  best  fit  to  the  data  by  an  ellipse,  the  Twiss  parameters 
were  determined.  The  values  of  the  beam  size  in  the  two  orthogonal  directions 
(X  and  Y)  and  the  respective  slopes  (.Y'  and  V  )  were  obtained  as:  .Y  =  14.5 
mm.  .Y'  =  132.2  mrad;  Y  =  18.9  mm.  Y'  =  175.4  mrad.  In  order  to  estimate 
the  beam  parameters  at  the  extraction  slit,  a  linear  beam  optics  code  is  used. 
This  code  takes  the  above  values  of  the  beam  parameters  at  r  =  11.75  cm  as 
input,  and  it  integrates  the  Kapchinskij-VIadimirskij  (K-V)  envelope  equations 
backwards  towards  the  extraction  aperture.  The  background  gas  may  have  an 
important  role  in  this  analysis,  particularly  in  relation  to  the  evaluation  of  the 
space-charge  neutralization  factor  fy  which  enters  in  the  calculation  of  the  space- 
charge  force.  X'arious  spatial  profiles  of  fy  have  been  considered:  the  results  are 
rather  insensitive  to  the  variation  of  fy  .  VVe  have  shown  two  cases  in  Fig.  2:  (a) 
fy  =  0,  and  (b)  fy  =  l.O  at  the  extraction  aperture  and  it  falls  off  exponentially 
afterwards.  From  this  analysis  we  predict  a  beam  radius  of  1.1  mm  and  slope  of 
72  mrad  at  the  extraction  slit. 

In  order  to  further  understand  the  role  of  a  neutralizing  background  on  mea¬ 
surements  of  the  spatial  evolution  of  beam  parameters,  we  studied  three  cases  of 
practical  relevance:  (a)  a  30  mA,  35  kV  beam  with  a  beam  radius  of  1.1  mm 
and  beam  divergence  of  72  mrad.  representing  a  magnetron  source  as  described 
above,  (b)  a  30  mA.  35  kV  beam  with  a  beam  radius  of  5.6  mm  and  near-parallel 
e.xtraction.  representing  a  BNL  volume  source.^  and  (c)  for  higher  current  -  120 
mA.  30  kV  beam  with  a  beam  radius  of  2  mm  and  near-parallel  extraction.  The 
effective  emittance  values  in  the  three  cases  are  almost  similar,  the  range  being 
5.0  — 6.6  X  10~^  m-rad.  Figure  3  shows  the  spatial  evolution  of  the  beam  envelopes 
when  several  different  values  of  the  space-charge  neutralization  factor  (fy  =  0. 
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Figure  1:  Phase-space  distribution  of  H~  beam  from  the  SSCL  magnetron  source. 
The  measurements  were  made  at  11.75  cm  downstream  from  the  e.xtraction  aper¬ 
ture. 

0.5,  1.0,  and  1.04)  were  assumed.  .As  a  typical  e.xampie.  fs  is  taken  constant,  C, 
up  to  z  =  10  cm.  and  then  it  linearly  drops  to  zero  at  z  =  15  cm.  The  results 
suggest  that  the  space-charge  neutralization  factor  significantly  changes  the  beam 
envelope  for  the  ctises  (b)  and  (c),  and  that  the  elfect  is  not  as  strong  in  the  case 
(a).  Comparing  the  various  terms  in  the  K-V  envelope  equations,  we  note  that 
the  ratio  of  the  space  charge  force  to  the  emittance  force.  Ka^/e^.  is  initially  (at 
z  =  0)  much  larger  than  1  in  the  cases  (b)  and  (c);  while  this  ratio  is  close  to 
1  in  the  case  (a).  ([\  =  lit,! is  the  generalized  beam  perveance.  A  is  the 
beam  current,  /q  is  the  characteristic  current  (3.1  x  10^  .A  for  H~  beam),  ^  =  u/c, 
7  =  (1  —  a  is  the  initial  beam  radius,  and  e  is  the  effective  beam  emit¬ 

tance.)  In  the  cases  where  the  space-charge  force  is  the  most  dominant  factor, 
it  may  be  postulated  that  the  aforementioned  method  to  estimate  initial  beam 
parameters  may  introduce  some  ambiguity,  unless  a  good  understanding  of  the 
gas  dynamics  is  achieved  and  is  incorporated  properly  in  the  analysis. 

SIMULATION  OF  BEAM  DYNAMICS  THROUGH  A  LEBT  SYSTEM 


(a)  H~  Beam  from  the  SSC  Magnetron  Source 

A  proto-type  ESQ  LEBT  system  has  been  developed  in-house  at  .Maryland, 
and  the  system  has  been  described  in  detail  elsewhere.^  Figure  4  shows  the  geom¬ 
etry  of  the  LEBT  system.  In  a  previous  article.^  we  reported  that  the  aperture 
of  the  lenses  in  the  proto-type  ESQ  LEBT  system  needs  to  be  increased  by  a  fac¬ 
tor  of  two  to  accommodate  the  large  excursion  of  the  SSC  magnetron  H“  beam. 
Further  details  of  the  beam  dynamics  are  presented  here,  where  the  ESQ  lens 
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Figure  2:  K-V  envelope  solution  for  H~  beam  propagating  over  the  drift  space  be¬ 
tween  the  extraction  aperture  of  the  magnetron  source  and  the  slit  of  the  emittance 
diagnostic.  Here  ;  =  0  is  considered  as  the  location  of  the  slit  of  the  emittance 
diagnostic;  the  beam  envelope  is  calculated  backward  toward  the  extraction  slit 
which  is  here  at  c  =  11.75  cm.  Space-charge  neutralization  factor,  /v.  is  taken 
as:  (a)  zero  (upper  figure),  and  (b)  exponentially  falling  with  the  peak  value  as 
1.0  at  2  =  11,75  cm  (lower  figure).  The  inset  in  the  lower  figure  shows  the  spatiaJ 
profile  of  f,v. 
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Figure  3:  Spatial  evolution  of  the  beam  envelope  for  three  different  cases:  (a) 
beam  current  I  =  30  mA,  beam  voltage  V  =  35  kV,  initial  beam  radius  a  =  l.l 
mm,  initial  beam  divergence  r'|r=a  =  "2  nuad;  (b)  Ik  and  K  are  same  zis  in  (a),  a 
=  5.6  mm.  r'|r=a  =  0:  (c)  /(,  =  120  mA.  Vi  =  30  kV.  a  =  2.0  mm.  r'\r=a  =  0.  The 
bottom  figure  shows  the  nature  of  the  assumed  space-charge  neutralization  factor 
/,v.  The  four  different  curves  in  each  figure  correspond  to  different  values  of  the 
constant  C.  Going  from  the  leftmost  side  of  each  figure,  solid  line:  C  =  0.  next 
dashed  line:  C  =  0.5.  dotted  line;  C  =  l.O,  next  dashed  line:  C  =  1.04. 
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Figure  4:  The  LEBT  system. 

parameters,  especially  the  voltages,  are  further  optimized  and  also  the  design  of 
the  end  section  of  the  LEBT  consisting  of  a  short  einzel  lens  is  further  developed. 

Figure  5  shows  the  most  optimized  envelope  solution  that  we  have  obtained 
until  now  after  a  number  of  iterations  of  our  numerical  scheme.’’  The  maximum 
beam  excursion  in  one  plane  (about  17  mm)  is  found  significantly  larger  than 
in  the  other  orthogonal  plane  (about  10  mm).  Initially  we  obtained  an  envelope 
solution  with  symmetric  beam  excursions  in  both  ,V  and  K-planes:  however,  the 
emittance  growth  was  higher  in  this  situation.  This  is  an  artifact  of  nonlinear 
beam  dynamics.  .\s  a  rule  of  thumb,  nonlinear  effects,  which  contribute  to  the 
distortion  of  the  beam  optics  and  thus  enhance  the  beam  emittance.  are  not 
important  if  the  beam  does  not  occupy  more  than  about  75%  of  the  quadrupole 
aperture.  Hence  it  is  clear  from  Fig.  5  that  the  emittance  growth  should  primarily 
occur  in  the  ,V-plane.  The  phase-space  distribution  of  the  particles,  obtained 
from  the  modified  P.\RMILA  code  (a  2|-D  code)  is  shown  in  Fig.  6;  the  results 
conform  to  the  physical  picture  revealed  in  Fig.  5.  Here  a  beam  current  of  about 
25  mA  is  transported.  The  remaining  part  of  the  30  mA  input  beam  current 
contributes  enormously  (more  than  300%)  to  the  emittance  growth.  This  group  of 
particles,  which  behaves  like  a  “‘halo’’  and  occupies  a  large  phase-space  boundary, 
is  eliminated  by  a  beam  scraper.  The  emittance  growth  for  about  25  mA  of  beam 
current  is  estimated  to  be  about  50%  assuming  a  K-V  type  distribution  of  the 
input  beam  (Fig.  6);  the  estimated  emittance  growth  is  about  60%  for  a  Gaussian 
input  beam  as  shown  in  Fig.  7. 

From  the  modified  P.\RMILA  results  shown  in  Fig.  6  it  is  estimated  that  the 
output  beam  parameters  at  6  mm  downstream  from  the  last  lens  of  the  ESQ  LEBT 
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Figure  5:  Envelope  solution  for  the  ESQ  LEBT  section  using  SSCL  magnetron 
beam  parameters.  .\  hard-edge  type  e.xternal  focusing  function  k(z),  as  shown  in 
the  bottom  figure,  is  assumed. 
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Figure  6;  .Modified  P.\RMIL.A  results  on  phcise-space  distribution  of  particles 
through  the  ESQ  LEBT  section  (magnetron  beam  case).  Left-side  figure:  input 
beam  distribution  (K-V  type):  right-side  figure:  output  beam. 
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Figure  7;  Modified  PARMILA  results  for  the  same  case  as  in  Fig.  6  using  a 
Gaussian  input  beam  distribution. 

are:  -Y  =  4.68  mm.  Y  =  4.38  mm.  .Y'  =  —49  mrad,  and  Y'  —  — 18  mrad.  This 
beam  is  given  as  input  to  a  short,  single  einzel  lens  section  cascaded  to  the  ESQ 
LEBT.  The  einzel  lens  is  composed  of  three  electrodes  as  shown  in  Fig.  8  (top 
figure)  -  the  first  and  the  third  electrodes  are  at  ground  potential,  and  the  center 
electrode  is  at  a  high  negative  potential  close  to  the  beam  voltage.  In  our  experi¬ 
mental  configuration,  the  two  grounded  electrodes  of  the  einzel  lens  are  made  as 
common  elements,  respectively,  with  the  last  shunt  piate  of  the  ESQ  LEBT  and 
the  front  wall  of  the  RFQ.  The  output  from  the  einzel  lens  is  determined  using 
the  SN0W-2D  code.  The  output  beam  distribution  in  Fig.  8  shows  the  results 
at  a  downstream  location  of  1.0  cm  from  the  front  of  the  RFQ’s  ground  wail:  an 
effective  ellipse  representing  the  output  beam  in  Fig.  6  is  given  as  input  here. 
The  beam  parameters  at  this  location  match  well  with  the  acceptance  ellipse  of 
the  RFQ.  However,  if  the  input  vane  of  the  RFQ  is  located  further  away,  the  lens 
parameters  need  to  be  adjusted.  It  may  be  noted  that  the  einzel  lens  section  does 
not  introduce  any  emittance  growth.  Hence  the  overall  emittance  growth  in  the 
entire  LEBT  channel  is  about  50%. 

(b)  H“  Beam  from  a  Volume  Source 

The  beam  parameters  of  BNL  volume  source^  are  used  in  this  case.  The  proto¬ 
type  ESQ  LEBT  system  developed  at  Maryland  is  well  suited  to  transport  such 
beams.  Figure  9  shows  the  K-V  envelope  solution;  we  have  aimed  here  to  obtain 
a  beam  convergence  of  about  —20  mrad  at  the  output  cis  the  final  focusing  will 
be  done  by  an  einzel  lens. 

The  phase-space  distribution  of  particles  in  the  input  and  output  beams  for 
the  ESQ  LEBT  section  is  shown  in  Fig.  10.  The  output  beam  is  almost  free  from 
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Figure  8:  SNOW-2D  results  for  the  einzel  lens  section  using  the  output  beam  from 
the  ESQ  LEBT  (magnetron  case).  Top  figure:  beam  envelope  through  the  lens. 
Middle  figure:  phase-space  plot  of  the  input  beam.  Bottom  figure:  phase-space 
plot  of  the  output  beam. 


Figure  9;  Envelope  solution  tor  the  ESQ  LEBT  section  using  volume  source  pa¬ 
rameters. 

distortions,  and  the  emittance  growth  is  insignificant.  The  SNOVV-2D  code  results 
showing  the  behavior  of  the  beam  particles  as  they  move  from  the  end  of  the  ESQ 
LEBT  section  and  are  being  focused  by  the  einzel  lens,  are  shown  in  Fig.  11.  The 
output  beam  distribution  is  taken  here  at  a  distance  of  1.0  cm  downstream  from 
the  front  end  of  the  third  electrode  (at  ground  potential)  in  the  einzel  lens:  the 
third  electrode  should  be  considered  as  the  front  wall  of  the  RFQ  in  reality.  The 
emittance  growth  in  the  entire  LEBT  channel  is  very  small.  <15%.  here. 

CONCLUSIONS 

One  of  the  major  points  in  this  article  is  the  simulation  analysis  to  estimate 
beam  parameters  at  the  e.xtraction  aperture  from  the  measured  data  at  a  certain 
distance  away,  typically  about  10  cm  downstream.  It  is  noted  that  an  uncertainty 
in  the  analysis  generally  occurs  due  to  the  lack  of  precise  knowledge  about  the 
neutralizing  gas  background.  The  importance  of  this  issue  depends  significantly  on 
the  initial  conditions  of  the  beam,  primarily  the  ratio  of  A'a^/e^.  For  Ka^fe^  ~  I 
close  to  the  extraction  aperture,  as  in  the  case  of  the  SSC  magnetron  beam,  the 
spatial  evolution  of  the  beam  envelope  over  a  distance  of  about  10  cm  is  quite 
insensitive  to  the  v'ariation  of  the  space-charge  neutralization  factor,  The 
beam  envelope  is  found  to  be  significantly  affected  by  the  variation  of  fs  for 
An*/ »  1:  examples  are  given  for  a  30  mA  volume  source  beam  and  a  situation 
with  large  beam  current  (120  mA). 

The  beam  dynamics  through  a  LEBT  system,  comprising  of  6- ESQ  lenses 
and  one  einzel  lens,  have  been  studied  in  detail  using  simulation  codes.  The 
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Figure  10:  Modified  PARMILA  results  on  phase-space  distribution  of  particles 
through  the  ESQ  LEBT  section  (volume  source  case).  Left-side  figure:  input 
beam  distribution  (K-V  type):  right-side  figure:  output  beam. 
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Figure  11:  SN0W-2D  results  for  the  einzel  lens  section  using  the  output  beam 
from  the  ESQ  LEBT  (volume  source  case).  Top  figure:  phase-space  plot  of  the 
input  beam.  Bottom  figure:  phase-space  plot  of  the  output  beam. 


analyses  have  been  done  using  two  widely  different  input  beam  parameters,  and 
the  LEBT  systems  have  been  designed  accordingly.  The  present  LEBT  system  is 
very  compact,  and  it  has  a  good  tunability  to  control  the  emittance  growth.  In 
the  case  of  a  30  mA.  3.5  kV  SSC  magnetron  beam  when  the  input  beam  radius 
is  taken  as  1.1  mm  and  the  beam  divergence  is  72  mrad.  the  beam  into  the  RFQ 
has  an  emitcance  growth  of  about  .50  to  60%.  Here  a  scraper  in  the  ESQ  LEBT 
system  is  used  to  deliver  25  mA  of  output  beam  current,  which  appears  to  be  the 
nominal  SSC  injector  beam  current.  This  part  of  the  beam  conforms  to  a  good 
core  geometry,  and  hence  it  suffers  less  due  to  aberrations.  NVhile  dealing  with 
a  larger  beam  (5.6  mm  in  radius  and  near-parallel  e.xtraction)  corresponding  to 
the  BNL  volume  source  case,  the  LEBT  system  is  shown  to  behave  very  nicely, 
and  the  full  beam  (30  mA)  can  be  matched  into  the  RFQ  with  practically  no 
emittance  growth. 

The  present  LEBT  system  has  several  strong  features.  The  transport  section 
of  about  30  cm  length  creates  a  good  buffer  space  between  the  ion  source  and 
the  RFQ.  It  thus  insures  a  clean  operation  of  the  RFQ.  Several  tuning  knobs 
allow  us  to  handle  some  variation  of  the  input  beam  characteristics  and  maintain 
an  acceptable  level  of  emittance  growth.  Anderson  et  al.'  reported  intriguing 
e.xperimental  results  with  similar  LEBT  systems  bcised  on  a  combination  of  ESQ 
lenses  and  one  einzel  lens  (referred  as  a  ring  lens  in  their  paper);  near-parallel  input 
beams  with  brightness  about  an  order  of  magnitude  lower  than  in  the  presently 
investigated  cases  were  considered  there. 

A  comment  is  made  in  the  conte.xt  of  the  present  e.xperimental  effort.  The 
current  LEBT  system  will  satisfy  the  RFQ’s  matching  condition  approximately 
I  cm  downstream  from  the  RFQ  aperture.  With  intense,  high- brightness  beams 
the  problem  of  matching  will  be  formidable  if  the  drift  space  between  the  end  of 
the  LEBT  and  the  RFQ  entrance  is  too  large;  the  emittance  growth  may  need  to 
be  sacrificed  enormously. 
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ESQ  FOCUSING  FOR  AN  INTENSE,  HIGH-BRIGHTNESS  H"  BEAM: 
EMITTANCE  GROWTH  AND  ITS  REMEDY 

S.  K.  Guharay,  C.  K.  Allen.  M.  Re'ser 
University  of  Maryland,  College  Park,  MD  20742  U.S.A. 

Abstract 

A  simple,  novel  electrostatic  quadrupole  (ESQ)  lens  system  has  been  developed  to  transport 
and  focus  an  intense,  high-brightness  (normalized  beam  brightness  ~  10“  A/m*  rad*  )  H~ 
beam.  The  physics  of  emittance  growth  in  the  ESQ  transport  system  is  studied  in  detail  by 
computer  simulation  of  beam  dynamics,  when  the  influence  of  various  factors,  e.g.,  input  beam 
parameters,  aberrations,  misalignments,  etc.  has  been  exanuned.  Possible  methods  to  control 
the  emittance  growth  have  been  suggested. 

I.  Introduction 

Although  there  has  been  a  significant  progress  in  the  physics  and  technology  of  high-energy 
accelerators,  many  new  areas  of  research  are  still  emerging  primarily  due  to  a  variety  of  needs  of 
the  new  accelerators  and  its  applications,  e.g.,  e.xploration  of  fundamental  particles,  space  defense 
by  antiballistic  neutral  particle  beams,  radio-active  waste  transmutation,  etc.  It  is  recognized  that 
a  very  good  quality  beam,  defined  by  its  emittance,  is  required  in  a  high-energy  accelerator.  This  is 
particularly  true  in  the  case  of  new  accelerators,  e.g.,  the  Superconducting  Super  Collider.*  where 
the  emittance  budget  through  the  accelerator  chain  is  very  tight  due  to  the  stringent  requirement  for 
high  luminosity  (lO'^^cm  ‘s~* )  of  the  colliding  beams.  In  this  pursuit  this  group  at  the  University 
of  Maryland  h.is  been  engaged  over  the  last  several  years  in  the  development  of  an  efficient  low- 
energy  beam  transport  (LEBT)  section.  This  unit  lies  between  an  ion  source  and  the  first  stage 
of  acceleration  of  the  beam  (e.g.,  a  radio- frequency  quadrupole  (RFQ)  accelerator)  and  plays  an 
important  role  in  accelerator  research. 

The  LEBT  acts  as  a  phase  transformer  of  the  beam  from  the  ion  source  so  that  a  normally 
diverging  beam  from  an  ion  source  can  be  focused  with  a  high  convergence  angle  (~  —100  mrad) 
Into  an  IIFQ.  Also,  the  LEBT  creates  a  butfer  space  between  the  ion  source  and  the  RFQ;  this 
reduces  the  possibility  of  any  contamination  of  the  RFQ  due  to  cesium  or  gas  load  from  the  ion 
source.  In  spite  of  an  e.xhaustive  body  of  literature  in  this  field,  the  state-of-art  is  not  satisfactory. 
So  far  the  major  experimental  efforts  on  LEBT  systems  in  accelerators  have  been  made  with 
jolenoidal  lenses*  or  magnetic  quadrupoles^  supplemented  by  gas  focusing.  The  design  of  such 
systems  has  been  mainly  empirical  and  control  parameters  are  not  known  apriori  with  sufficient 
accuracy.  Significant  emittance  growth  of  the  beam  is  noted  in  this  scheme.  Other  schemes  deal  with 
electrostatic  lenses  (einzel  lens.'  electrostatic  quadrupole  lens.^®  helical  electrostatic  quadrupole 
lens' )  or  radio- frequency  quadrupole  lens.^  .Although  such  alternate  schemes  have  been  tested  in 
some  experiments,  no  e.xperimental  effort  has  been  made  so  far  in  regard  to  transporting  and 
focusing  intense,  high-brightness  H~  beams  (typically,  beam  current  of  .lO  mA,  beam  voltage  of  .15 
kV,  normalized  brightness  'v  io“  .A/m^rad^).  It  is  important  to  focus  our  attention  to  such  beams 
<as  they  are  relevant  to  the  potential  applications  of  modern  accelerators. 

The  present  article  delineates  the  problem  of  H"  beam  transport  and  its  matching  to  an  RFQ 
using  a  novel  6-lens  electrostatic  quadrupole  (ESQ)  system  as  an  LLBT.  The  main  emphasis  here  is 
to  understand  the  physics  of  emittance  growth  in  an  ESQ  transport  system  and  the  possible  ways 
to  mitigate  the  problem.  In  the  past  the  ESQ  LEBT  system  was  used  very  successfully  to  transport 
H  beams;®  however,  the  beam  perveance  was  much  lower  there.  Here  we  consider  a  Iiighly  space- 
charge  dominated,  high-brightness  beam  (beam  perveance  K  =  0.003  and  normalized  brightness 


--  iO“  A/m-rad*).  It  is  quite  a  challenging  task  to  transport  sucn  beams  over  a  certain  aistance  and 
focus  it  into  an  RFQ  without  any  significant  emlttance  dilution.  The  design  is  based  on  detailed 
computer  simulation  of  beam  dynamics  through  the  transport  section.  The  computer  predictions 
regarding  performance  of  the  LEBT  system  are  analyzed  in  the  context  of  beam  parameters  of  two 
particular  types  of  H~  ion  sources  -  Penning- Dudnikov  source  as  used  in  the  BE.AR  experiment.^ 
and  a  volume  ionization  source  of  Brookhaven  National  Laboratory.^  These  two  types  of  sources 
are  particularly  suitable  to  obtain  high- brightness  H”  beams. In  our  previous  articles.  H”  beams 
from  a  Pennlng-Dudnikov  source  were  mainly  considered.  Here,  the  emphasis  is  on  the  transport 
of  H~  beams  from  a  volume  source.  We  have  also  addressed  the  issue  of  sensitivity  of  the  beam 
parameters  due  to  misalignments.  Finally,  the  problem  of  matching  the  H~  beam  to  an  RFQ  is 
discussed. 

II.  Beam  Dynamics  through  the  ESQ  LEBT 

The  design  procedure  as  reported  earlier  **  involves  a  sequence  of  computer  code  analysis:  (i) 
.A  linear  beam  optics  code  integrates  the  K-V  envelope  equations  and  generates  the  basic  geo¬ 
metrical  parameters  of  the  lens  system,  (ii)  .\  3-D  LAPL.ACE  solver  calculates  the  equipotentials 
and  evaluates  the  fringe-field  matrices  as  suggested  by  Matsuda  and  WoUnik.'^  (iii)  .A  modified 
PAR.MIL.A  code^^  uses  input  from  the  above  two  steps  and  evaluates  the  beam  parameters  (beam 
size,  emittance  growth,  etc.).  This  scheme  led  us  to  choose  a  combination  of  six  ESQ  lens  system 
to  transport  H~  beams  over  a  distance  of  30  cm  and  to  provide  a  moderate  convergence  (~  -  >0 
mrad)  at  the  end.  The  technical  details  of  th^  LEBT  svstem  and  some  results  on  ilie  computer 
code  predictions  of  beam  c.iaracteriitics  have  been  given  earlier."  The  essential  points  to  note  from 
our  previous  articles  are: 

1.  .\  unique  feature  in  the  design  of  the  ESQ  LEBT  system  is  that  the  entire  system  is  seif- 
aligned  mechanically,  riiis  demands  high  precision  in  the  fabrication  of  the  components  and 
their  assembly.  On  the  other  hand,  it  eliminates  ilie  usually  tedious  alignment  job  in  the 
'•xperiment. 

1.  The  performance  of  the  ESQ  LEBT  has  been  examined  using  beam  parameters  of  a  Penning- 
Dudnikov  source.  Some  assumptions  of  the  input  beam  (normalized  brightness  =  x  10‘“ 
A/m*rad*.  divergence  at  full  beam  radius  of  I  ;>’m  =  JO  mrad)  are  made  on  the  basis  oi 
analysis  of  the  emittance  data  at  -  =  10.6  cm  from  the  extraction  slit."  The  output  beam 
from  the  ESQ  LEBT  is  found  converging  (~  -25  mrad).  and  it  shows  an  emittance  growth  by 
a  factor  of  about  l.S  assuming  a  I\-V'  type  distribution  of  the  beam.  The  oiriittance  growtii 
is  identified  as  due  mainly  to  hromatic  aberrations. 

The  above  analysis  has  been  carried  out  further.  Figure  1(a)  shows  the  evolution  of  rms  nor¬ 
malized  emittance  through  the  ESQ  transport  channel.  I'he  enhancement  of  emittance  occurs 
essentially  in  the  second  and  fifth  lenses,  when  the  amplitude  of  the  beam  envelope  grows  to  more 
than  ''0%  of  the  ((uadrupole  aperture.  The  beam  particles  residing  in  the  outer  jiart  of  the  enve¬ 
lope  are  responsible  for  the  emittance  growth,  and  these  particles  can  be  rejected  by  suitable  use 
of  the  ground  plates  between  the  adjacent  ESQ  lenses  (Fig.  10  in  ref.  1 1 )  as  beam  scrapers.  Tao 
beam  scrapers  are  inserted,  one  in  front  of  the  second  lens  and  the  other  in  front  of  the  fifth  lens. 
Figure  Kb)  shows  the  evolution  of  emittance  through  the  ESQ  LEBT  when  10%  of  beam  particles 
are  scraped  out.  The  emittance  growth  is  a  factor  of  l.l  here.  This  is  an  efficient  ’.  ay  to  deliver 


an  almost  emittance-preserved  beam  to  an  RFQ.  The  above  analysis  has  been  carried  out  using 
typical  parameters  of  H“  beams  from  a  volume  ionization  source,  here  a  BNL-type.  Usually  the 
beam  from  such  sources  has  lower  current  density  compared  to,  say,  the  Penning-Dudnikov  case. 
Hence,  a  larger  e.xtraction  slit  (radius  of  extr?ction  slit  =  5.6  mm  here  as  opposed  to  a  radius  of 
1  mm  In  the  Penning-Oudnikov  case)  is  used  to  draw  the  same  amount  of  current  in  the  case  of 
volume  sources.  Also  the  beam  divergence  at  the  extraction  slit  can  be  very  small:  almost  a  parallel 
beam  car.  be  extracted  from  the  source.  The  estimated  normalized  brightness  of  this  beam  is  about 
3  X  10‘°.\/ra*rad^;  all  other  parameters  have  been  given  earlier.'^ 

Figure  2  shows  the  modified  PARMILA  results  on  distribution  of  the  beam  particles  correspond¬ 
ing  to  an  H~  beam  from  the  volume  source.  Here  the  parameters  corresponding  to  a  4  times  rms 
ellipse,  constructed  by  including  90%  of  the  beam  particles  at  the  output  of  the  ESQ  LEBT,  are: 
■Vm**  =  3.3  mm.  =  2.7  mm,  (A'a,*,)'  =  -13.5  mrad,  {YmmxY  =  -14.9  mrad.  The  maximum 
voltage  on  the  ESQ  lenses  was  4.7  kV.  The  emittance  growth  of  the  beam  is  negligibly  small,  ~ 
1%.  Tliis  result  suggests  that  the  characteristics  of  the  input  beam  are  very  crucial  to  the  issue  of 
•■mittance  growth  in  a  transport  system  -  -  a  parallel  input  beam  is  desirable. 

The  sensitivity  of  output  beam  parameters  with  variation  of  beam  voltage,  beam  current,  mis- 
uienment  of  the  beam  axis  with  respect  to  the  LEBT  system  has  been  studied.  The  K-V  code 
Kiaivsis  .suggests  that  the  beam  parameters  do  not  change  noticeably  for  a  ±1%  change  of  quad 
voltage  (from  the  ideal  setpoint)  on  all  the  lenses  simultaneously.  Similar  insensitivity  is  noted 
lor  variation  of  beam  current  (ideal-  =  -30  mA)  within  a  few  milliamperes.  Preliminary  analysis  on 
the  aforementioned  misalignments  is  done  using  the  modified  P.\RMILA  code  (no  image  charge  is 
included  in  the  present  code).  With  the  variation  of  the  amount  of  off-centering  of  the  input  beam, 
the  phase-space  distribution  of  the  output  beam  remains  almost  invariant,  while  the  beam  centroid 
-hifts  coherently,  .V  translational  off-centering  by  1  mil  at  the  input  is  amplified  by  about  a  factor 
of  2  off-contering  at  the  output  of  the  LEBT:  this  also  introduces  an  angular  error  of  the  beam 
•  oMtroid  by  about  2  mrad.  An  on-centered  input  beam  with  an  error  of  1  mrad  in  the  injection 
ancle  shows  an  off-centering  of  the  output  beam  by  about  l.S  mil:  the  error  in  the  angle  siavs 
•dmosi  the  same  through  the  LEBT  channel. 

III.  Matching  to  an  RFQ 

I  he  beam  trom  the  LEBT  section  needs  to  be  matched  with  the  RFQ  in  order  to  achieve  a  good 
transmission  and  preserve  emittance  through  this  channel.  The  typical  Twiss  parameters  for  the 
acceptance  ellipse  of  an  RFQ  (e.g..  DEAR  and  SSC  RFQ)  demand  a  beam  convergence  of  about  -90 
mrad  at  the  full  beam  radius  of  about  1.3  mm.  It  is  an  arduous  task  to  satisfy  these  requirements 
of  an  RFQ  by  an  ESQ  LEBT  alone  without  sacrificing  the  emittance  growth.  .An  LEBT  composed 
of  two  modules  —  the  ESQ  lenses  as  a  beam  transport  section  and  an  einzel  lens  at  the  end  as  a 
matching  section  —  seems  to  be  a  good  choice  in  this  pursuit. 

figure  .)  shows  a  schematic  of  the  ESQ  LEBT  with  a  short  modular  attachment  containing 
an  einzel  lens.  Preliminary  analysis  is  done  using  the  K-V  code.  The  beam  envelopes  for  the 
two  situations  are  shown  in  Fig. 4:  (a)  Penning- Dudnikov  case,  (b)  volume  source  case.  The  beam 
parameters  at  the  end  of  the  einzel  lens  match  well  with  the  requirements  of  the  RFQ  as  mentioned 
earlier.  .A  high  value  (about  7.0)  of  the  ratio  of  the  aperture  of  the  einzel  lens  to  the  beam  size  is 
taken  here:  hence,  the  einzel  lens  is  expected  to  have  an  insignificant  contribution  to  the  emittance 
growth.  This  problem  is  being  currently  investigated  to  get  a  quantitative  answer. 


rV.  Experiments  and  Dbcussion 

The  ESQ  LEBT  system  has  been  constructed  in-house.  The  overall  mechanical  alignment  of  the 
apparatus  has  been  measured,  and  found  to  be  within  ±  1.5  mil.  The  power  supply  system  for  the 
ESQ  lenses  has  been  attached.  Detailed  voltage  hold-off  tests  have  oeen  done,  and  the  performance 
is  satisfactory.  It  is  planned  to  test  this  ESQ  LEBT  with  an  H~  beam  from  a  magnetron  source 
at  the  SSC  Laboratory.  The  code  results  are  being  reviewed  in  the  light  of  the  beam  parameters 
from  the  magnetron  source. 

It  has  been  shown  that  a  combination  of  an  ESQ  LEBT  system  with  an  einzel  lens  can  be  used 
very  effectively  to  transport  a  space-charge  dominated,  high- brightness  H“  beam  over  a  distance 
of  about  30  cm  and  focus  the  beam  into  an  RFQ  without  any  significant  emittance  dilution.  This 
method  appears  to  be  attractive  particularly  to  handle  a  highly  diverging  beam  from  an  ion  source. 
Such  a  scheme  has  a  number  of  advantages:  (i)  flexibility  to  handle  a  wide  range  of  input  beam 
parameters  from  different  types  of  ion  sources,  (ii)  allow  sufficient  buffer  space  (not  field-free) 
between  an  ion  source  and  an  RFQ  for  differential  pumping,  (iii)  experiments  with  low-voltage 
power  supplies  over  the  major  part  of  the  transport,  and  (iv)  ease  of  fine  tuning.  It  is  now  a  very 
important  issue  to  determine  the  optimum  length  of  the  I.EBT  section  to  insure  a  reliable  operation 
of  the  injector  system  in  an  accelerator.  Carefully  controlled  experiments  with  various  transport 
schemes  are  warranted  and  a  good  database  is  required  to  unravel  the  reliability  ol  the  computer 
simulation  results. 
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Fig.  1.  Evolution  of  rms  normalized  emittance  (solid  line:  X-component;  dashed  line: 
Y-component)  through  the  ESQ  lenses  for  an  H“  beam  from  the  Penning- Dudnikov 
source:  (a)  without  any  beam  scraper;  fb)  with  beam  scrapers,  one  in  front  of  the 
second  lens  and  the  other  in  front  of  the  fifth  lens.  The  numbers  inside  the  figure 
identify  the  location  of  the  downstream  endpoint  of  the  corresponding  lens. 


Fig.  2.  .Modified 
PARMILA  results  of  distri¬ 
bution  of  beam  particles  for 
an  H”  beam  from  the  BNL 
volume  source:  Input  beam 
(top  figure!:  Output  beam 
(bottom  figure). 


I  i^.  4.  K-\  envelope  .solution  (solid  lino:  .K  eomponent;  dashed  line:  \’-(Oinponeni ) 
for  the  matched  beam  to  an  RFQ  using  Jf“  beam  parameters  corresponding  to:  (a) 
I  enning- Dudni kov  source:  (b)  BNf.  volume  source. 
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ABSTRACT 

A  fast  approximation  method  to  the  3D  electrostatic  problem  is 
developed.  The  method  of  moments  procedure  is  outlined  for  the  particular 
case  of  Laplace's  equation.  The  resulting  matrix-vector  equation  is  solved  by 
a  conjugate  gradient  algorithm.  These  techniques  are  then  implemented  with 
a  computer  code  running  on  a  PC  and  used  to  solve  example  problems. 

1.  INTRODUCTION 

We  are  interested  in  the  transport  of  low  energy  ions  beams, 
specifically,  we  have  been  engaged  in  the  design  and  analysis  of  Low  Energy 
Beam  Transport  (LEBT)  section  of  an  accelerator  column’  *.  To  avoid  gas 
focusing  we  use  electrostatic  lenses  for  the  relatively  slow  moving  ions.  The 
action  of  such  lenses  is  completely  characterized  by  their  spatial  potential 
distribution  <{>.  Thus,  analysis  of  electrostatic  lenses  invariably  requires  the 
solution  of  Laplace's  equation. 

Although  simple  in  form,  this  equation  must  be  solved  numerically  for 
most  geometries  of  practical  interest.  Many  approximation  techniques  are  very 
successful  to  this  end  (e.g.  finite  differences)  and  are  covered  extensively  in  the 
literature.  However,  for  a  fully  3D  treatment  computation  and  machine  storage 
usually  become  extreme.  Manipulation  of  the  solution  data  also  becomes  quite 
cumbersome.  Consequently,  these  situations  require  special  hardware  such 
as  a  supercomputer.  We  present  an  approximation  technique  which  is  fully  3D 
and  remarkably  computational  efficient.  Ideally,  we  wish  to  implement  this 
technique  as  CADware  for  the  IBM  PC. 

The  technique  relies  on  a  combination  of  the  method  of  moments  and 
fast  iterative  techniques  for  solving  linear  systems.  Specifically,  we  proceed  by 
reformulating  Laplace's  equation  into  an  integral  equation  over  the  boundary 
surface,  reducing  the  dimensionality  of  the  original  system.  The  new  problem 
is  approximated  by  the  method  of  moments*  to  yield  the  matrix-vector  equation 
Ax=y.  We  then  use  conjugate  gradient  algorithms^  to  solve  this  equation. 


’Supported  by  ONR/SDIO 


2.0  PROBLEM  FORMULATION 


In  this  section  we  clarify  the 
physical  system  which  we  wish  to  model. 

Figure  1  depicts  an  abstract  geometric 
representation  of  the  electrostatic 
problem.  We  have  a  closed  surface  F  in 
Euclidean  3-space  which  represents  the 
boundary  of  our  problem.  F  separates  E^ 
into  two  regions  Q,  and  Q.  representing 
the  (bounded)  interior  of  F  and  the 
(unbounded)  exterior  of  F.  respectively.  Figure  1:  Electrostatic  Problem 
Typically  we  associate  the  union  of  F  and 

Qj  as  a  conductor  in  3-space  while  we  are  interested  in  the  potential  4»  in  the 
region  external  to  the  conductor.  This  is  an  example  of  the  so  called 
exterior  Dirichlet  problem  and  can  be  classified  mathematically  as  follows: 

V®<l)(x)  =  0  VxeQ^, 

<|)(x)  =  VxeF. 

The  function  f  is  the  given  set  of  boundary  values  and  it  constitutes  the  data  of 
the  problem.  For  a  conductor  f  is  a  constant  over  the  boundary  F.  In  addition, 
we  usually  impose  the  physical  restriction  that  (t>  vanishes  at  infinity. 

Also  of  interest  is  the  interior  Dirichlet  problem  where  we  are  concerned 
with  the  distribution  (j)  in  the  region  Q,.  This  situation  is  formulated 
mathematically  below. 


V^(t)(x)  =0  VxeO, ,  ^2) 

4)(x)  =  /(x)  VxeF. 

Notice  that  both  problems  take  the  same  boundary  values,  namely  f. '  Typically 
for  the  interior  problem  we  also  require  that  (j)  remains  bounded  in  n^. 

We  may  reformulate  both  of  these  problems  into  a  single  integral 
equation  on  F*.  The  alternate  problem  is  achieved  by  the  introduction  of  an 
intermediate  function  o  defined  on  F. 


(5) 


<t)(x)  =  jG(x,x)a(x')dx^  VX€Q 
r 

where  Q  is  the  union  of  Qj  and  Q,.  Thus  equation  (5)  is  valid  in  all  space 
except  the  boundary  r.  We  recognize  G(x,x')  as  the  free  space  Green's 
function  for  Poisson's  equation  and  equation  (5)  as  the  potential  due  to  a 
surface  charge  density  a  in  free  space.  It  is  interesting  to  note  that  even 
though  problems  (1)  and  (2)  are  not  related  in  any  obvious  manner,  solution  of 
(3)  yields  the  correct  o  in  both  cases  given  the  same  boundary  values  f. 

3.0  THE  METHOD  OF  MOMENTS 

We  now  approximate  o  in  equation  (3)  using  the  method  of  moments. 
It  will  be  convenient  to  rewrite  (3)  by  defining  the  operator  K.  We  have 

Ka  AfG(x,x')a(x')dx' =  f  fe/.p(r)  (6) 

r 

where  we  have  chosen  f  in  Up  in  ordej-  to  analyze  convergence.  We  see  that 
K  is  a  compact  linear  self-adjoint  integral  operator  with  kernel  G(x,x').  Since  F 
is  compact,  we  may  choose  a  countable  set  of  functions  {u„}  which  are  dense 
in  Lp(n,  that  is  they  form  a  basis  for  the  vector  space  Lp(r)-  These  functions 
are  usually  referred  to  as  the  set  of  expansion  functions.  Accordingly,  any 
element  of  Lp  may  be  approximated  arbitrarily  well  by  c  inear  combination  of 
these  functions.  Thus, 


for  some  set  {aj  of  coefficients.  By  substituting  (7)  into  (6)  and  using  the 
linearity  of  K  we  have 

E  a^u„  =  f.  (8) 

/f-1 

Now  we  choose  another  set  of  functions  {vj  in  L,  (where  1/p  +  1/q  =  1)  to  be 
used  as  a  set  of  weighting  functions.  By  taking  the  inner  product  of  (8)  with 
each  of  the  v„,'s  we  get  a  set  of  equations  of  the  form 

m 

^  a„iKu„,vJ  =  {f,v„)  m=1,2,3,...  (9) 

ff-i 

where  <  -,  -  )  denotes  the  usual  inner  product 

-  {u,v)  =  Ju(x)  v(x)0&f.  (iO) 

r 


Note  that  since  f  is  given  and  {u„}  and  {v^}  are  selected  arbitrarily,  the  only 
unknowns  are  the  coefficients  {aj.  To  obtain  an  approximation  for  o  we  must 
truncate  the  indices  n  and  m  at  a  finite  value.  If  we  choose  the  same  value  for 
both  indices,  say  N,  we  are  then  left  with  the  following  matrix-vector  equation. 

:  :  =  I  (11) 

which  can  be  written  more  compactly  as 

/lx  =  y.  (12) 

The  matrix  A  is  known  in  the  literature  as  the  moment  matrix  for  (6),  since  we 
take  the  moments  of  Ka  with  respect  to  the  weights  {vj.  Equation  (12)  may 
be  solved  by  standard  matrix  techniques  to  yield  a  solution  for  x,  the  vector  of 
coefficients  (  a,,  ...  ,  a^  We  then  have  our  approximation  to  o  according  to 
equation  (7). 

COMMENTS 

(1 )  The  method  of  moments  generates  solutions  that  converge  in  the  mean,  i.e. 
in  the  Lp  norm.  The  exact  Lp  space  where  convergence  occurs  depends  on 
the  choice  of  expansion  and  weighting  functions. 

(2)  When  we  chose  {vj  =  {u  J,  we  have  Galerkin's  method  which  is  known  to 
be  equivalent  to  the  Rayleigh-Ritz  variational  method^  Thus,  we  see  that 
the  method  of  moments  is  a  generalization  of  the  Rayleigh-Ritz  procedure. 
For  Galerkin's  method  we  have  convergence  in  the  norm  (mean-squared 
convergence),  where  the  approximate  a  lies  in  span  {u„}. 

(3)  By  choosing  {v  J  =  {5(x-x„,)}  where  6  is  the  Dirac  delta  function  and  the  x,„'s 
are  some  set  of  points  on  r,  we  have  a  point-matched  solution.  Such 
solutions  are  known  to  converge  in  the  L,  norm  (pointwise  convergence). 

(4)  For  valid  solutions  we  must  make  sure  the  set  {KUp}  spans  the  range  of  K. 
Otherwise,  we  converge  to  solutions  to  the  problem  Ko  =  Pf  where  P  is  the 
projection  operator  onto  the  space  span{Ku„}. 

(5)  K  ’  is  unbounded  and  there  exists  f  in  Lp  such  that  (8)  has  no  solution  o  in 
Lp  (such  solution  may  however  be  interpreted  distributional'y).  However,  K 
is  a  positive  operator,  indeed  <Ka,o>  is  recognized  as  the  electrostatic 
energy  in  the  system.  Consequently,  in  the  discrete  approximation,  0  is  not 
an  eigenvalue  of  A  (it  is  in  the  limit  N->a),  though). 

4.0  CONJUGATE  GRADIENT  ALGORITHM 

We  now  turn  our  attention  to  the  task  of  solving  the  matrix-vector 
equation  (12).  The  traditional  approach  is  to  use  some  direct  method  such  as 


LU  decomposition  or  Gaussian  elimination.  With  these  techniques  the  amount 
of  computation  necessary  to  solve  an  N"*  order  system  is  known  a  priori.  For 
example,  it  is  known  that  inversion  of  (12)  requires  O(N^)  operations.  Instead, 
we  have  chosen  an  iterative  technique  where  the  amount  of  computation  is  not 
known  in  advance. 

The  conjugate  gradient  algorithms  are  a  specialization  of  the  more 
general  technique  of  conjugate  directions  methods.  These  techniques  are 
expressly  developed  to  solve  the  problem® 

min  =  ^x^Ax  -  y^x.  (13) 

It  is  assumed  that  the  matrix  A  is  positive  definite  so  that  a  solution  does  exist. 
Note  that  the  solution  to  this  problem,  obtained  by  setting  the  gradient  of  g  to 
zero,  is  given  by  x  =  A'V-  Conjugate  direction  methods  are  based  on  the  idea 
of  generating  a  complete  set  of  linearly  independent  vectors  {d„}  which  have 
the  property  d/Ad„,  =  0  whenever  m?tn.  This  is  known  as  A-orthogonality,  or 
A-conjugacy. 

Instead  of  solving  (12)  directly,  in  conjugate  gradients  we  choose  to 
iteratively  minimize  some  functional,  for  example  g  in  the  above  equation.  The 
value  of  X  furnishing  this  minimum  is  the  solution  to  (12).  The  method  starts 
with  an  initial  guess  for  x^,  then  it  generates  a  sequence  X|  that  minimizes  the 
functional.  The  sequence  {x}  will  converge  to  the  exact  solution  x  in  a  finite 
number  of  iterations.  Even  if  the  matrix  A  is  not  invertible  in  the  classical  sense 
(0  is  an  eigenvalue  of  A),  the  conjugate  gradient  algorithm  will  converge  to  a 
solution  in  the  least  squares  sense.  The  algorithm  we  have  chosen  was  taken 
from  Sarkar  el.  ai*  The  functional  which  it  minimizes  is  the  (Z^)  norm  squared 
of  the  residual  r,  =  y-Axj  (i.e.  minimize  ||y-AXjf  ). 

5.0  IMPLEMENTATION 


To  apply  the  method  of  moments  it  is 
necessary  that  we  select  an  appropriate  set 
of  expansion  functions  {u„}  and  weighting 
functions  {v„,}.  Since  our  aim  is  to  minimize 
computation,  we  have  selected  a  point¬ 
matching  procedure.  Thus,  our  weighting 
functions  {v„}  are  given  by  the  set  {5(x-xJ} 
where  the  x„'s  are  the  match  points  in  r  to 
be  determined.  This  technique  yields  the 
least  computing  time,  and  the  procedure  is 
straightforward.  It  also  yields  good  results 
as  long  as  one  is  careful  in  the  selection  of 
the  match  point  locations^*. 

We  choose  a  set  of  piecewise 
constant  functions  for  {u„}.  First,  f  is 


Figure  2;  Triangulated  Sphere 


divided  into  triangular  subdomains,  denoted  by  the  set  {TJ  where  n  runs  from 
1  to  N.  If  r  has  curvature  then  we  have  an  approximation  to  r  as  well  (for 
example  figure  2  depicts  the  triangulated  approximation  to  a  sphere).  Now  let 
the  u„'s  be  a  set  of  piecewise  constant  functions  on  each  T„. 


1  if  xeT„ 
0 


(14) 


It  can  be  shown  that  such  a  set  is  dense  in  Lp  so  a  may  be  approximated 
arbitrarily  well®.  Thus,  according  to  comment  (4)  we  must  be  certain  that  Pf  = 
{f(x„)}  is  in  the  span  of  {Ku„}.  We  shall  see  later  that  this  criterion  places 
restrictions  on  our  triangulation  of  F.  We  have  chosen  the  matching  point  set 
{x„}  to  be  the  centroids  of  each  triangle  T„.  Since  this  is  the  center  of  mass  for 
our  constant  expansion  functions,  it  seems  to  be  a  reasonable  approximation 
to  f  across  the  triangle  face. 

In  order  to  form  equation  (11)  we  must  determine  the  inner  products 
<Ku„,5(x-xJ>.  For  our  choice  of  expansion  functions,  this  amounts  to  the 
evaluation  of  the  following  integral. 


=  /  G{x„x)dx' 


_1_ 

4tc 


_ 1 _ 

/{x„-xf*(y„-yf^{z„-zf 


dx’dy'dz' 


(15) 


This  can  be  done  completely 
numerically,  hybrid  analytic  and 
numeric,  or  completely 
analytically’®. 

6.0  EXAMPLES 

We  have  implemented  the 
above  techniques  in  a  computer 
program  written  in  Borland  C++. 

The  platform  is  an  i486  PC 
operating  at  33  MHz  and  running 
Windows  3.1  operating  system. 

All  examples  were  run  in  single 
precision  arithmetic  except  where 
noted.  We  have  taken  arbitrary  units. 


I  I - nuTiwc  . arial-rc 

Figure  3;  Conducting  Sphere  Axial 
Potential 


CONDUCTING  SPHERE 

The  analytic  potential  distribution  for  a  conduction  sphere  of  radius  1  and  held 


at  a  potential  1  is  known  to  be 


1  if  \x\< 

=  1 

-L  //  |xl  <:  1. 

i  lx! 

Figure  3  shows  a  comparison  of  this  analytic  formula  with  the  numerical  results 
achieved  for  a  triangulated  sphere  of  432  triangles  (shown  in  figure  2).  The 
solution  converges  to  a  norm  squared  residual  error  <  10‘®  in  75  iterations 
taking  about  7  minutes  real  time  (double  precision  required  42  iterations  within 
5  minutes). 

A  single  number  which  is  indicative  of  solution  quality  is  the  capacitance.  This 
value  may  be  calculated  by  recognizing  o  as  surface  charge  density.  We 
numerically  calculated  the  capacitance  to  be  110.2  pF  for  the  above  situation 
in  MKS  units,  the  true  value  is  111.3  pF. 


EINZEL  LENS 

An  einzel  lens  was  modeled  as  two 
cylindrical  pipes  of  radius  5  and 
length  10,  separated  axially  by  a 
distance  of  3.  Both  pipes  were 
capped  with  plates  having  an  inner 
radius  of  1/2.  A  conducting  plate 
with  outer  radius  5  and  inner  radius 
1/2  was  centered  between  the  two 
pipes  at  position  z=0.  The  center 
plate  was  driven  to  a  potential  of 
1  while  the  two  outer  pipes  were 
grounded.  EI-Kareh  provides  an 
approximate  analytic  axial  potential 
distribution  given  by” 

4,(2)  (2z+3)  ’  (2z-3)  -2zfar7*’(2z)].  (17) 

3n  2  2 

This  approximation  assumes  three  infinite  plates  of  inner  radius  1/2  with 
uniform  fields  at  2  =  ±00.  The  discretization  was  made  with  464  triangles.  This 
system  took  261  iterations  to  converge  to  the  usual  error  criterion  (about  18 
minutes  of  real  time).  What  is  happening  here  is  that  our  boundary  function  f 
of  equation  (8)  is  almost  outside  the  range  of  {Ku„}.  Fortunately,  our  iterative 
technique  is  indicating  this  situation  by  slower  convergence!  In  order  to  avoid 
such  situations  it  is  necessary  to  reevaluate  the  triangulation  of  the  system  to 
insure  that  {Ku„}  spans  the  space  of  boundary  functions.  Even  though  our 
patch  model  was  borderline  ill-conditioned,  the  stability  of  the  conjugate 


Figure  4:  Einzel  Lens  Axial  Potential 


gradient  method  still  yielded  reasonable  results  (see  figure  4). 


7.0  CONCLUSION 

The  method  outlined  in  this  paper  exhibits  enough  computational 
efficiency  to  allow  practical  implementation  on  a  PC.  It  yields  accurate  results 
in  a  reasonable  amount  of  real  computing  time.  Moreover,  the  conjugate 
gradient  algorithm  furnishes  suitable  error  criterion  to  judge  the  quality  of  the 
solution.  However,  the  method  does  require  more  work  to  implement  than,  say, 
the  rather  straightforward  technique  of  finite  differencing. 
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Abstract 

A  low-energy  beam  transport  system  is  designed 
with  the  aim  of  transporting  a  30  mA.  35  kV  H" 
beam  from  a  volume  source  and  focusing  it  into  an 
RFQ.  The  characteristics  of  the  beam  from  the  source 
are  determined  analyzing  the  emittance  data.  The 
behavior  of  the  beam  through  the  LEBT  is  stud¬ 
ied  using  simulation  codes.  The  system  parameters 
are  optimized  so  that  the  LEBT  has  a  very  modest 
contribution  to  the  emittance  growth  (here  a  factor 
of  about  1.5)  and  the  emittance  budget  of  the  linac 
section  is  maintained. 

I.  INTRODUCTION 

One  of  the  vital  considerations  in  modern  high- 
energy  accelerators  is  related  to  the  design  of  an  effi¬ 
cient  low-energy  beam  transport  (LEBT)  section  so 
that  an  intense,  high-brightness  beam  (here  we  con¬ 
sider  an  H“  beam)  can  be  transported  over  certain 
distance  and  finally  focused  into  the  commonly  used 
RFQ  accelerator  in  the  linac  section.  The  emittance 
growth  in  the  LEBT  is  the  key  issue  in  developing 
a  good  scheme  at  the  low-energy  end  of  the  acceler¬ 
ator  chain.  Thus,  in  order  to  achieve  a  good  beam 
quality  and  match  it  to  the  acceptance  of  an  RFQ  a 
systematic  study  of  beam  dynamics  in  the  preceed- 
ing  sections  including  the  e.xtraction  optics  of  the 
ion  source  is  warranted.  VVe  have  e.xperienced  that 
emittance  measurements  of  the  beam  from  an  ion 
source  and  an  analysis  of  the  data  to  characterize 
the  beam  at  the  extraction  electrode  are  two  impor¬ 
tant  problems  in  the  context  of  designing  an  efficient 
LEBT.  Earlier  we  reported  (1,2)  on  the  beam  char¬ 
acterization  and  beam  dynamics  through  a  LEBT 
for  H“  beams  from  a  Penning-Dudnikov  type  (3) 
and  magnetron  type  sources.  In  recent  years  sig¬ 
nificant  progress  has  been  made  in  the  performance 
of  volume  sources  (4,5);  normalized  beam  brightness 
approaches  about  10"  A/(m-rad)*  for  H“  beams. 
This  article  highlights  on  the  study  of  H~  beams 
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Figure  1:  Contour  emittance  plot. 

from  the  SSC  volume  source  with  the  aim  of  design¬ 
ing  a  LEBT  system  to  deliver  a  30  mA.  35  kV  beam 
matched  with  the  RFQ  input. 

II.  BEAM  CHARACTERISTICS 

A.  H~  Beam  from  the  Ion  Source 

The  H“  beam  from  the  SSCL  volume  source  is 
measured  at  10.13  cm  downstream  after  the  electrons 
(ratio  of  initial  electron  to  ion  current  40)  are  de¬ 
flected  away  from  the  extracted  beam  current  by  a  10 
cm  long  magnetic  trap.  Figure  1  shows  the  contour 
plots  in  the  z  -  x'  space  from  emittance  diagnos¬ 
tics;  the  flattening  of  the  distribution  in  the  upper 
half  is  possibly  caused  by  the  space-charge  force  due 
to  the  electrons  deflected  upward.  The  beam  pa¬ 
rameters  at  X  =  10.13  cm  are:  beam  size  D  =  2.38 
cm.  full  divergence  Ad  =  260  mrad.  ;rf„  =  0.1537a- 
mm-mrad.  These  data  are  used  in  an  envelope  sim¬ 
ulation  code  to  estimate  the  beam  characteristics  at 
the  extraction  electrode;  space-charge  effects  due  to 
the  electrons  in  the  extraction  region  are  included. 
Figure  2(a)  shows  the  assumed  space-charge  correc¬ 
tion  factor.  /,  due  to  the  electrons.  Note  that  /  is 
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Figure  '2;  (a)  Space-charge  correction  factor  f;  (b) 
beam  envelope,  z  ~  10.13  cm  corresponds  to  the 
location  of  emittance  measurements. 

negative  here  and  the  beam  perveance  is  to  be  mul¬ 
tiplied  by  the  factor  (I  -  /).  The  beam  envelope  in 
Fig.  J(b)  is  evaluated  by  integrating  the  K-V  en- 
vplope  equations  using  a  fourth-order  Runge-Kutta 
method.  This  analysis  suggests  that  the  beam  at  the 
'■xtraction  electrode  emerges  nearly  parallel,  and  the 
beam  size  is  close  to  the  aperture  radius  (=  4  mm). 

B.  Desired  Output  Beam  Parameters  from  the  ESQ 
LEBT 

The  purpose  of  the  LEBT  section  is  to  isolate 
the  RFQ  from  the  ion  source  for  a  clean  operation 
and  also  to  deliver  a  matched  beam  to  the  RFQ.  The 
SSC  RFQ  acceptance  for  a  30  mA,  35  kV  H“  beam 
is  given  by  the  Twiss  parameters:  a  =  1.26,  0  = 
1.86  cm/rad.  =  0.2  t  mm-mrad.  As  the  normal¬ 
ized  rms  emittance  of  the  H~  beam  from  the  source 
is  about  0.1537  ir  mm-mrad  (Fig.  1).  the  LEBT  is 
to  be  designed  under  a  very  tight  emittance  bud¬ 
get.  The  matching  condition  dictates  that  the  beam 
parameters  at  the  tip  of  the  RFQ  vane  should  be: 
beam  radius  =  1.3  mm  and  the  corresponding  slope 
of  the  beam  envelope  =  -89  mrad;  this  is  located  at 
about  3  cm  downstream  from  the  front  wall  of  the 


RFQ.  It  has  been  shown  earlier  that  a  short  (about  5 
cm  long)  sineie  einzel  lens  module  between  an  ESQ 
LEBT  and  the  RFQ  will  be  a  good  choice  in  satis¬ 
fying  the  aforementioned  stringent  conditions  of  the 
RFQ  (2].  The  ESQ  LEBT  transforms  the  highly  di¬ 
verging  beam  from  the  ion  source  into  a  moderately 
converging  one  without  any  significant  emittance  di¬ 
lution.  and  the  einzel  lens  provides  the  final  strong 
focusing.  This  analysis  showed  that  the  parameters 
of  the  output  beam  from  the  ESQ  LEBT  should  fol¬ 
low:  beam  radius  3  to  5  mm.  the  corresponding 
slope  of  the  beam  envelope  ~  -30  to  -50  mrad. 

III.  BEAM  TRANSPORT  THROUGH  THE  ESQ 
LEBT 

The  design  principles  of  the  ESQ  LEBT  follow 
the  scheme  as  discussed  earlier  (1.2].  The  present 
configuration  of  the  magnetic  trap  in  the  extraction 
region  of  the  SSCL  volume  source  restricts  the  ESQ 
LEBT’s  distance  to  the  extraction  aperture  to  about 
10  cm.  This  causes  the  beam  to  blow  up  significantly 
(Fig.l).  After  a  detailed  analysis  with  such  a  beam 
it  is  recognized  that  the  goal  to  deliver  a  matched 
beam  to  the  RFQ  for  the  full  beam  current  (30  mA) 
is  a  very  difficult  task.  Our  analysis  suggests  that 
a  shorter  magnetic  trap  (about  5  cm  long)  will  be 
a  better  choice.  Figure  3  shows  the  beam  envelope 
through  the  ESQ  LEBT  when  a  hard-edge  focusing 
function  for  the  external  field  is  assumed.  An  initial 
drift  space  of  5  cm  long  is  considered,  and  a  pro¬ 
file  of  the  space-charge  correction  factor  due  to  the 
electrons  (Fig.  3,  bottom)  is  assumed.  The  beam  pa¬ 
rameters  at  the  extraction  aperture  are  taken  from 
the  analysis  of  Fig.  2.  The  maximum  aperture  ra¬ 
dius  of  the  quadrupoles  is  22  mm;  it  was  taken  as 
12  mm  when  the  LEBT  was  closer  (1.5  cm)  to  the 
extractor  (ij. 

The  distribution  of  the  beam  particles  through 
the  ESQ  LEBT  is  estimated  using  a  modified  P.\R- 
MILA  code  (6).  Figure  4  shows  the  particle  dis¬ 
tribution  in  phase  space  for  /  =  -30  mA.  The  esti¬ 
mated  output  beam  parameters  are:  .V  =  3.5  mm. 
V'  =  3.5  mm.  .V'  =  -51.7  mrad.  Y'  =  -50.8  mrad. 

The  emittance  growth  is  primarily  due 
to  chromatic  aberrations. 

IV.  CONCLUSIONS 

Emittance  measurements  of  a  30  mA.  35  kV  H~ 
beam  from  the  SSCL  volume  source  have  been  stud¬ 
ied.  The  simulation  results  suggest  that  the  H~ 


f.h-  ■.micij, 


o*- 


;  ccmi 

Figure  3:  K-V  envelope  solution. 


beam  envelope  has  a  waist  at  the  extraction  aper¬ 
ture.  With  this  definition  ol  the  input  beam  and 
a  given  set  of  characteristic  parameters  of  the  RFQ 
acceptance,  we  have  designed  a  LEBT  system.  The 
particular  design  of  the  LEBT  consisting  of  six  ESQ 
lenses  and  a  short  einzei  lens  can  transport  the  full 
beam  current  and  match  it  to  the  RFQ. 

Beam  transport  experiments  with  a  prototvpe 
ESQ  LEBT  will  be  conducted  to  validate  the  simu¬ 
lation  predictions.  Further,  the  3-D  LAPLACE  sim¬ 
ulation  scheme  is  being  improved  using  a  method  of 
moments  where  any  arbitrary  boundarv  can  be  rep¬ 
resented  numerically  and  the  practical  problems  will 
be  simulated  more  realistically  [7]. 

This  work  was  supported  bv  ONR/SDIO  and 
DOE. 
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Abstract 

A  compact.  6-lens  electrostatic  quadrupole  (ESQ) 
LEBT  system  has  been  constructed  at  the  University 
of  Maryland  to  transport  a  30  mA.  35  kV  H"  beam 
over  a  distance  of  about  30  cm.  A  short  einzei  lens 
section  is  included  at  the  end  of  the  ESQ  LEBT  to 
establish  a  good  matching  of  the  beam  to  the  radio¬ 
frequency  quadrupole  (RFQ)  accelerator,  and  to  meet 
the  emittance  requirements  of  the  linac  in  the  Super¬ 
conducting  Super  Collider.  Computer  code  predic¬ 
tions  on  the  beam  dynamics  through  the  LEBT  with 
cxperimentaily  measured  input  beam  data  are  dis-" 
cussed. 

Introduction 

An  elficient  ion  source-cum-low  energy  beam  trans¬ 
port  I  LEBT)  section  is  highly  desired  to  deliver  a 
good-quality  beam  to  the  low-energy  booster  (LEB) 
in  the  collider  ring  chain  of  SSCL.  The  intrinsic  emit¬ 
tance  of  the  ll~  beam  from  an  ion  source,  volume- 
ionization  or  magnetron  type,  is  typically  about  0.12 
T  mm-mrad  ( rms  normalized  value);  the  LEB  requires 
that  the  transverse  beam  emittance  at  the  output  of 
the  linac  section  be  <  0.3jr  mm-mrad  (!].  The  compo¬ 
nents  of  the  linac  between  the  ion  source  and  the  LEB 
.ire:  LEBT.  radio- frequency  quadrupole  (RFQ)  ac¬ 
celerator.  drift-tube  linac  (DTL)  and  coupled-cavity 
linac  ( CCL).  Computer  code  analyses  of  the  beam  dy¬ 
namics  through  the  OTL  and  the  CCL  suggest  that 
the  emittance  growth  in  these  two  sections  is  not  sig¬ 
nificant.  being  in  the  range  of  10-15%  [2].  The  per¬ 
formance  of  the  RFQ,  e.g.,  beam  transmission  and 
emittance  growth,  depends  primarily  on  the  Twiss 
parameters  of  the  beam  at  its  input:  the  transverse 
rms  normalized  emittance  of  the  input  beam  is  de¬ 
sired  to  be  <  0.2jr  mm-mrad.  Hence  the  LEBT’s 
contribution  to  the  emittance  budget  must  be  main¬ 
tained  within  a  factor  of  1.6  of  the  input  beam  emit¬ 
tance.  This  suggests  that  the  ion  source-cum-LEBT 
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Figure  1:  The  LEBT  system. 


section  plays  a  critical  role  in  the  performance  of  the 
linac. 

This  article  addresses  some  important  problems  in 
designing  an  efficient  LEBT  system.  Here,  an  ESQ 
lens  system  is  primarily  considered:  two  other  vari¬ 
ants  of  the  LEBT  system,  einzei  lenses  and  a  helical 
quadrupole  lens,  are  also  being  investigated  in  the 
present  context  at  SSCL  [.3].  The  analyses  are  based 
on  computer  code  simulations,  where  the  input  beam 
parameters  are  mostly  taken  from  experimental  mea¬ 
surements. 

Beam  Dynamics,  Design  of  LEBT  and 
Discussions 

In  our  previous  paper  [4],  we  described  a  (J-lens 
ESQ  LEBT  system  developed  .it  the  L'niversity  of 
.Maryland.  The  predicted  performance  of  the  ESQ 
LEBT  are  now  examined  in  the  light  of  beam  param¬ 
eters  relevant  to  the  SSCL  program.  Two  types  of  H“ 
sources  are  considered  in  the  SSCL  injector  develop¬ 
ment  -  a  volume  source  and  a  magnetron  source. 

In  the  context  of  the  volume  source,  we  have  used 
the  beam  parameters  corresponding  to  a  Brookhaven 
.Vational  Laboratory  {BNL)-type  source  (5].  .A 
30  mA.  35  kV  H"  beam  is  extracted  through  a  I  cm^ 
circular  aperture.  .A  parallel  beam  is  assumed  to 
emerge  from  the  extraction  aperture.  The  beam  en¬ 
velope  through  a  compact  6-lens  ESQ  LEBT  section 
(Fig.l),  which  has  been  constructed  in-house  at 


Figure  2:  Linear  beam  optics  result  for  a  parallel  in¬ 
put  beam.  Top:  Amplitude  (X,  Y)  of  the  beam  enve¬ 
lope  through  the  ESQ  LEBT  along  .x  (solid  line!  and 
y  (dotted  line);  bottom:  focusing  function,  k.  r  is  the 
direction  of  propagation  of  the  beam. 

Maryland,  is  computed  by  integrating  the  K-V  en¬ 
velope  equations  using  a  fourth-order  llunge-Kutta 
method.  Figure  2  shows  the  beam  envelopes,  when  a 
linear  external  focusing  force  represented  by  a  hard- 
edge  focusing  function.  n{:),  is  assumed.  .\3  men¬ 
tioned  in  the  previous  article  (dj,  the  ii.SCL  RFQ  re¬ 
quires  a  circular  beam  of  about  l.J  mm  in  radius  and 
a  beam  convergence  of  about  -90  mrad  at  To 

match  these  conditions  without  sacrificing  the  emit- 
tance  growth,  an  additional  unit,  a  single-piece  einzel 
lens  module,  is  included  at  the  end  of  the  ESQ  LEBT 
section  (Fig.l).  The  particle  distribution  at  the  out¬ 
put  of  the  ESQ  LEBT  is  computed  using  a  modified 
P.\RMILA  code  [6],  and  it  does  not  show  any  signif¬ 
icant  emittance  growth,  ~  ')%.  Results  in  a  similar 
situation  have  been  shown  earlier  (•(]:  this  point  is  not 
elaborated  further.  The  einzel  lens  turns  the  moder¬ 
ately  convergent  (~  —20  mrad)  beam  from  the  ESQ 
LEBT  into  a  strongly  convergent  (~  -110  mrad) 
beam  with  a  negligible  emittance  dilution.  ~  -5%.  The 
behavior  of  the  beam  through  the  einzel  lens  section, 
predicted  by  the  SNOW-2D  code,  is  shown  in  Fig.3. 

The  aforementioned  analyses  have  been  carried  out 
using  the  parameters  of  a  30  mA.  35  kV  11“  beam  e.\- 
tracted  from  the  SSCL  magnetron  source.  The  char¬ 
acteristics  of  the  beam  are  measured  at  a  distance  of 
11.75  cm  downstream  from  the  tip  of  the  extraction 
cone.  Using  these  results  as  initial  beam  conditions, 
the  beam  parameters  at  the  tip  of  the  extraction  cone 
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Figure  3;  .S.XO\V-2D  results  for  the  einzel  lens  mod¬ 
ule.  lop:  Beam  trajectory  through  the  einzel  lens 
module.  The  center  electrode  is  at  -35  kV  and  the 
two  end  electrodes  are  grounded.  .Middle:  Phase- 
space  distribution  of  the  input  beam.  Bottom:  Phase- 
space  distribution  of  the  output  l)eain. 

arc  predicted.  .\  plausible  estimate  ol  the  beam  pa¬ 
rameters  at  the  tip  of  the  extr.aciion  cone  is:  beam 
radius  =1.1  mm.  and  divergence  at  Tma*  =  "2  mrad. 

The  lens  aperture  in  the  ESQ  LEBT  in  Fig.l  is  not 
large  enough  to  accommoilate  the  highly  diverging 
11“  beam  from  the  magnetron  source.  .\  preliminary 
design  study  reveals  that  the  aperture  of  the  KSQ 
lenses,  second  through  fifth  in  Fig.l.  i.s  to  l>e  inrre.a.sed 
by  a  factor  of  2:  this  demands  a  higher  voltage  on  the 
quadrupoies.  The  beam  dynamics  through  the  KSQ 
LEBT  is  followed  using  the  moilifieji  I’.MIMIL.V  code. 
Figure  4  shows  the  output  beam  distribution  y  bottom 
figure),  when  a  K-V  type  input  beam  (lop  figure)  is 
.assumed.  The  ground  plate  in  front  of  the  second  lens 
in  Fig.l  has  been  used  .as  a  beam  scraper  to  reject 
.about  15%  of  the  beam  particles,  which  contribute 
significantly  to  the  emittance  growth.  The  output 
beam  in  Fig.4  still  suffers  from  some  distortions,  giv¬ 
ing  rise  to  an  emittance  growth  by  a  factor  of  about 
1.5.  Further  optimization  of  the  ESQ  LEBT  is  war¬ 
ranted  to  improve  the  present  situation.  .Nevertheless. 
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Figure  4:  Modified  PARMILA  results  on  particle  dis¬ 
tribution  for  H“  beam  from  the  SSCL  magnetron 
source.  Top:  input  to  the  ESQ  LEBT.  Bottom:  out¬ 
put  from  the  ESQ  LEBT. 

an  output  beam  current  of  ~  25  m.A.,  as  required  for 
the  SSC  RFQ,  is  achievable  from  the  ESQ  LEBT.  In 
regard  to  matching  the  beam  to  the  RFQ,  an  einzel 
lens  module  is  included  at  the  end  of  the  ESQ  LEBT. 
Figure  5  shows  preliminary  results  of  beam  transport 
through  the  einzel  lens  module.  The  smooth  nature 
of  the  input  beam  is  an  artifact  of  modeling  it  from 
an  estimate  of  the  effective  values  of  beam  parame-' 
ters  in  Fig.  4.  The  einzel  lens  does  not  contribute 
to  the  emit.ance  growth:  the  beam  parameters  at  the 
front  end  of  the  third  electrode  (at  ground  potential) 
match  closely  to  the  acceptance  ellipse  of  the  SSCL 
IIFQ. 

The  above  analyses  for  the  LEBT  have  been  done  in 
two  separate  stages-  first  the  ESQ  lenses  with  modi¬ 
fied  P.\R.\1ILA  and  then,  the  einzel  lens  with  S.NOW- 
JD.  Tlie  simulation  predictions  will  hold  well,  if  the 
matching  between  the  two  codes  is  complimentary. 
An  effort  is  being  made  to  include  the  option  of  an 
••inzel  lens  in  the  modified  P.\RMIL.\.  when  a  more 
reliable  design  tool  will  be  available. 

Conclusions 

riie  problem  of  low-energy  beam  transport  aiul  its 
matching  to  an  RFQ  has  been  studied  in  reference  to 
two  special  cases  of  the  input  H~  beam:  (i)  a  parallel 
beam,  and  (ii)  a  highly  divergent  beam.  The  compu¬ 
tational  results  suggest  that  the  key  point  in  design¬ 
ing  an  efficient  H~  injector  (ion  source-cum-LEBT) 
relates  to  achieving  a  well-conditioned  beam,  e.g.,  un¬ 
aberrated  and  near- parallel,  from  the  ion  source.  .An 
emphasis  is  latid  here  to  develope  an  LEBT  appara¬ 
tus  with  compactness,  mechanical  stability  and  fle.x- 
ibility  for  easy  modifications.  .\  combination  of  a 
6-lens  ESQ  module  and  one  short  einzel  lens  module, 
as  adapted  in  the  present  LEBT  design,  appears  to 
be  a  good  choice  in  this  respect.  An  experiment  is 
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Figure  5:  S.N'0\V-2D  results  for  the  einzel  lens  mod¬ 
ule.  Top:  Beam  trajectory.  The  center  olcctrouc  is 
at  -36  kV.  .Middle:  Phase-space  distribution  of  the 
input  beam  ( Effective  values  from  Fig.  I  arc  used  to 
model  it.)  Bottom:  Phase-sp.icc  distribution  of  the 
output  beam. 

bring  planned  on  the  SSCL  test  stand  to  stinlv  the 
beam  characteristics  through  the  ESQ  LEBT  system 
developed  at  .Maryland,  and  test  the  reliability  of  sim¬ 
ulation  predictions. 
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Abstract 

With  an  aim  of  transporting  an  initially  diverging  high-perveance  (generalized  beam  per- 
veance  2/(,//o/3^7^  =  0.003),  high- brightness  (normalized  brightness  ~  10“  A/(m-rad)*)  H~ 
beam  and  finally  focusing  it  without  any  significant  emittance  dilution,  a  detailed  simulation 

scheme  has  been  set  up  incorporating  the  various  nonlinear  forces  due  to  the  beam  and  the 

• 

external  focusing  elements,  e.g.,  due  to  space  charges,  geometrical  and  chromatic  aberrations. 
The  analysis  is  done  following  a  particular  hierarchy  to  identify  the  mechanism  of  emittance 
growth;  this  procedure  is  used  to  optimize  the  lens  parameters.  A  combination  of  six  elec¬ 
trostatic  quadrupole  lenses  is  configured  to  deliver  a  satisfactory  solution.  The  estimated 
emittance  growth  is  a  factor  of  about  1.6,  and  this  is  mainly  due  to  chromatic  aberrations.  A 
relatively  small  group  of  particles  is  found  to  be  responsible  for  the  emittance  growth.  The 
analysis  highlights  a  number  of  important  issues,  e.g.,  sensitivity  to  the  beam  distribution, 
beam  current,  lens  misalignments,  etc.  An  ESQ  LEBT  system  with  some  novel  features  in 
terms  of  compactness  and  mechanical  rigidity  is  developed,  and  its  essential  characteristics 
are  described. 
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1.  Introduction 


y. 

The  study  of  high- brightness  charged  particle  beam  transport  has  great  relevance  in 
many  modern  applications.  In  today’s  and  next  generation’s  high  energy  colliders,  e.g., 
Tevatron,  SSC,  NLC,  etc.,  one  of  the  vital  requirements  is  to  achieve  luminosity  of  colliding 
charged  particle  beams  of  order  10^^  cm~^  s~'  or  higher;  this  demands  that  beam  brightness 
also  be  very  high  [1].  H"  be2uns,  with  normalized  brightness  of  ~  10'*  A/(m-rad)*,  are 
required  in  space  defense  for  generation  of  intense  particle  beams  to  probe  any  foreign  objects. 
The  importance  of  high  intensity,  high-brightness  beams  is  adso  evident  in  heavy- ion  fusion 
(HIF),  free  electron  lasers,  etc.,  and  of  late,  in  an  attractive  scheme  for  radio-active  waste 
transmutation  [2].  The  recent  trends  of  ion-beam  related  research  in  accelerators  and  fusion 
reveal  that  major  activities  have  been  initiated  with  H“  beams  due  to  their  merits  over 
proton  beams: 

(a)  Given  a  phase-space  area,  the  intensity  of  an  H~  beam  can  be  enhanced  significantly 
by  adapting  the  principle  of  non-Liouvillian  stacking.  This  is  particularly  useful  in 
modern  accelerators. 

(b)  The  charge-exchange  cross  section  at  energy  ~  100  keV  is  much  higher  for  H~  beams; 
hence,  H“  beams  are  used  for  the  development  of  high-energy  neutral  beams  in  mag¬ 
netic  fusion  research  and  also  in  space  defense. 

During  its  long  travel  through  the  various  components  in  an  accelerator,  the  character¬ 
istics  of  an  ion  beam  are  signiflcantly  determined  by  its  behavior  in  the  early  stages,  e.g., 
extraction  optics  in  the  ion  source  and  subsequently,  the  low-energy  beam  transport  (LEBT) 
section.  The  crux  of  the  problem  is  to  obtain  a  high  quality  beam  (usually  defined  by  the 
beam  emittance)  from  an  appropriate  ion  source  with  minimum  distortions  due  to  aberra¬ 
tions  and  nonlinear  forces,  and  to  preserve  its  quality,  as  much  2is  possible,  in  the  transport 
and  acceleration  chain.  This  article  focuses  on  the  LEBT  section  as  a  systematic  study  on 
this  part  is  still  missing,  particularly  in  relevance  to  transport  and  focusing  of  space-charge 
dominated,  high-brightness  H~  beams. 
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The  LEBT  section  isolates  the  first  stage  of  acceleration,  e.g.,  a  radio-frequency  c^adrupole 
accelerator  (RFQ),  from  an  ion  source  providing  a  buffer  space  for  differential  pumping;  also, 
it  presents  a  clean,  matched  beam  to  an  RFQ.  The  gas  focusing  scheme,  usually  supplemented 
by  magnetic  quadrupoles  or  solenoidal  magnets,  has  been  mostly  used  as  a  LEBT  system 
[3-5].  The  advantage  of  the  scheme  is  that  the  experimental  hardware  is  simple,  and  it  is 
capable  of  handling  large  beam  current,  ~  100  mA.  However,  this  scheme  is  not  well  suited 
to  beam  pulses  shorter  than  the  gas  neutralization  time,  typically  >  50  ns.  Also,  due  to  the 
complex  nature  of  atomic  and  molecular  processes  involved,  it  is  difficult  to  set  up  accurate 
computer  simulation  models  and  obtain  reliable  predictions  apriori.  Thus  our  understanding 
of  this  approach  remains  at  a  qualitative  level,  and  there  is  some  ambiguity  in  the  definition 
of  the  control  parameters  in  experiments.  One  of  the  most  impressive  applications  of  the  gas 
focusing  scheme  has  been  in  the  context  of  the  BEAR  (Beam  Experiment  Aboard  a  Rocket) 
accelerator  program  at  Los  Alamos  [4].  Results  on  the  sensitivity  of  beam  parameters  to 
variations  of  the  gas  pressure  suggest  that  a  reproducible  operation  of  such  a  scheme  is  very 
difficult.  The  other  LEBT  schemes  include  radio-frequency  quadrupole  (RFQ)  lenses  [6]  and 
electrostatic  lenses  [7-12],  e.g,,  einzel  lens,  electrostatic  quadrupole  (ESQ)  lens,  and  helical 
electrostatic  quadrupole  (HESQ)  lens.  Until  now,  the  utility  of  RFQ  lenses  hcis  been  least 
explored.  In  the  category  of  electrostatic  lenses,  einzel  lenses  are  commonly  used.  The  con¬ 
ventional  einzel  lenses  have  a  drawback  in  that  plasmas  may  build  up  in  the  field-free  region, 
rendering  the  system  susceptable  to  beam-plasma  instabilities.  Anderson’s  “ring”  lens  ver¬ 
sion  appears  to  be  free  from  this  problem  [8|.  The  einzel  lenses  require  power  supplies  close 
to  the  beam  voltage,  typically  ~  30  kV  for  a  35  kV  beam;  the  ESQ  lenses  have  an  edge  over 
the  einzel  lenses  in  this  respect.  Several  practical  factors,  e.g.,  tuning  with  low-voltage  power 
supplies  (typically  <10  kV),  elimination  of  any  field- free  region,  and  ease  of  computer  mod¬ 
eling,  make  the  ESQ  system  a  very  attractive  choice  as  a  LEBT.  Reiser  [13]  reported  briefly 
a  comparative  study  of  the  aforementioned  LEBT  schemes.  It  is  understood  that  in  spite  of 
a  good  deal  of  research  involving  LEBT  systems,  a  systematic  study  of  beam  dynamics  and 
optimization  of  LEBT  parameters  has  not  been  made  in  the  context  of  transporting  intense, 
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high-brightness  H~  beams.  In  order  to  obtain  insight  into  the  problem  and  advance  the 
present  state-of-the  art  of  LEBT  systems,  detailed  simulation  studies  should  be  first  made, 
and  the  simulation  predictions  compared  with  experimental  results.  In  a  recent  article  Bru 
enumerated  a  beam  transport  program  in  which  ail  the  forces  are  considered  linear  [14].  The 
present  article  addresses  in  detail  for  the  first  time  the  beam  dynamics  issues  relevant  to 
designing  zm  ESQ  LEBT  for  intense,  high-brightness  H“  beams  incorporating  the  various 
nonlinear  effects  due  to  aberrations,  fringe-fields,  etc.;  a  special  emphasis  is  given  here  to 
understanding  emittance  growth  and  methods  to  control  it.  Typical  beam  parameters  of  a 
Penning-Dudnikov  type  ion  source  are  considered  in  which  the  generalized  beam  perveance 
K  =  0.003  and  normalized  beam  brightness  Bn  ~  10^'  A/(m-rad)^).  Here,  A'  is  defined  as 
the  ratio  of  2/i/ where  A  is  the  beam  current,  Iq  is  the  charaicteristic  current  (3.1  x  10^ 
A  for  H”  beam),  0  =  vjc,  7  =  (1  — 

Some  key  physics  issues  on  emittance  growth  in  an  ESQ  LEBT  and  its  control  are  in¬ 
vestigated  in  detail.  The  simulation  scheme  is  developed  in  steps.  A  linear  beam  optics 
code  is  written  to  integrate  the  well-known  K-V  envelope  equations.  The  lens  parameters 
are  approximated  from  this  analysis.  A  3D  Laplace  solver  maps  the  equipotentials  of  the 
lens  system,  and  the  fringe  fields  are  evaluated.  Several  articles  have  dealt  with  fringe  fields 
[15-20].  Following  the  method  of  Matsuda  and  Wollnik  [19],  the  effect  of  the  fringe  fields  is 
included  here  in  a  particle  simulation  code,  which  essentially  is  a  modified  version  of  the  well- 
known  PARMILA  code  [21].  The  evolution  of  emittance  through  the  ESQ  LEBT  channel  is 
determined  from  the  modified  PARMILA  code.  This  simulation  scheme  is  iterated  until  a 
satisfactory  solution  is  obtained.  The  analysis  is  made  here  in  the  context  of  transporting 
an  initially  round  (radius  =  I  nun),  diverging  (slope  =  20  mrad  at  the  full  beaun  radius)  30 
mA,  35  kV  H~  beam  over  a  length  of  about  30  cm  and  transforming  it  into  a  round  (radius 
~  1  mm),  converging  (slope  ~  —40  mrad)  beam.  The  control  parameters  of  an  ESQ  lens 
can  be  simulated  reaisonably  well,  and  this  permits  to  delineate  an  in-depth  investigation  of 
the  various  sources  of  emittance  dilution.  This  study  thus  establishes  a  strong  foundation 
for  designing  an  efficient  LEBT  system.  Some  comments  are  made  on  the  development  of 
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the  LEBT  system  and  its  experimental  tests;  details  are  beyond  the  scope  of  this  paper  and 
these  issues  will  be  reported  elsewhere. 

Section  2  pertains  to  the  beam  dynamics.  Some  characteristic  features  of  the  LEBT 
system  are  given  in  Section  3.  Section  4  includes  conclusions. 

2.  Beam  Dynamics  through  the  ESQ  LEBT 
A.  Linear  beam  optics  calculations 

The  study  of  beam  dynamics  in  an  ESQ  LEBT  system  is  developed  from  a  rather  simple, 
idealistic  model  in  order  to  set  up  the  lens  parameters  grossly.  Afterwards,  the  simulation  is 
done  including  various  practical  features  ,  e.g.,  real  geometry  of  the  electrodes,  fringe  fields, 
aberrations,  etc.  First,  a  beam  optics  code,  which  includes  a  linear  external  focusing  force, 
here  a  hard-edge  type,  in  the  K-V  envelope  equations  is  used;  this  integrates  the  following 
coupled  equations  by  fourth-order  Runge-Kutta  method  and  gives  the  behavior  of  the  beam 
envelope. 

2K  c* 

X''  +  K,{z)X-  3^  =  0, 

'>/<■  C? 

-  yty  " 

Cj.  and  ty  are  the  x  and  y  components  of  the  unnormalized  beam  emittance,  respectively. 
The  focusing  function  due  to  the  externally  applied  force  is  k  =  V^/{Vi,R^),  where  K,  is  the 
applied  voltage  on  the  quadrupole,  H  is  the  beam  extraction  voltage,  and  is  the  aperture 
radius  of  the  aperture  radius  of  the  quadrupole.  The  derivatives  are  taken  with  respect 
to  z,  which  follows  the  direction  of  propagation  of  the  beam.  The  amplitude  of  the  beam 
envelope,  X(z)  and  Y{z),  is  determined  using  certain  values  of  the  input  beam  parameters 
and  an  assumed  configuration  of  the  lens  geometry.  In  the  case  of  any  gas  accumulation 
in  the  channel,  which  may  occur  at  the  interface  between  an  ion  source  and  a  LEBT,  the 
above  equations  are  solved  with  the  beam  perveance  term  K  multiplied  by  (1  —  7^/);  /  is 
the  charge  neutralization  factor  and  7  =  (I  — 

In  order  to  validate  the  results  predicted  by  the  linear  beam  optics  code,  the  various 
nonlinear  contributions,  namely  due  to  aberrations,  fringe  fields  and  image  effects,  are  to  be 
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minimized.  The  following  constraints  control  the  nonlinear  effects  {22|. 

(i)  If  the  maximum  beam  excursion  is  not  allowed  to  exceed  10%  of  the  length  of  a 
quadrupole  /,  spherical  aberrations  may  be  kept  at  low  level. 

(ii)  To  reduce  chromatic  aberrations,  the  change  in  the  beam  energy  due  to  the  quadrupole 
focusing  field  should  be  less  than  about  5%  of  the  total  beam  energy 

(iii)  The  image  field  may  be  neglected  if  the  maximum  excursion  of  the  beam  envelope 
remains  within  75%  of  the  quadrupole  aperture  radius  R^. 

The  following  empirical  relationship  is  satisfied  to  avoid  any  voltage  breakdown  [27] 

d  (cm)  >  1.4  X  (kV), 

where  d  is  the  interelectrode  spacing  and  V  is  the  voltage  difference. 

Due  to  limited  accessibility  in  the  neighborhood  of  the  beam  extraction  region,  it  is 
difficult  to  obtain  a  reliable  knowledge  of  beam  parameters  at  the  extraction  point;  this 
introduces  some  ambiguity  in  the  initial  values  of  the  beam  parameters  which  are  given  as 
input  to  the  code.  Usually  measurements  are  made  at  a  distance  of  about  10  cm  downstream 
from  the  extraction  slit  and  the  data  are  interpolated  using  some  simulation  model  to  predict 
beam  parameters  at  the  extraction  point.  Table  1  lists  the  estimated  values  of  the  H“  beam 
parameters  corresponding  to  the  Penning-Dudnikov  source  used  in  the  BEAR  experiment, 
when  the  nominal  emittance  data  [24]  measured  at  a  distance  of  10.6  cm  downstream  from 
the  extraction  aperture  have  been  used. 

The  desired  solution  is  obtained  after  several  iterations,  when  the  lens  parameters  (e.g., 
aperture  and  dimensions  of  the  individual  lenses,  spacing  between  the  lenses,  and  number 
of  lens  elements)  are  adjusted  at  each  step.  As  mentioned  earlier,  the  beam  envelope  is 
constrained  at  each  stage  through  the  LEBT  system  so  that  distortions  due  to  aberrations 
and  image  effects  are  minimized.  The  LEBT  system  has  been  thus  configured  using  a  set 
of  six  ESQ  lenses.  The  geometrical  parameters  of  the  lenses  are  given  in  Table  2.  Here  I  is 
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the  length  of  an  electrode  and  L  is  the  separation  between  adjacent  lenses.  Note  that  the 
geometrical  parameters  conform  to  a  symmetric  triplet  configuration. 

Figure  1  shows  the  K-V  envelope  solution  (top  figure),  when  the  input  beam  parameters 
correspond  to  Table  1  and  a  hard-edge  type  focusing  function  (bottom  figure),  k(2),  is 
assumed.  The  voltages  on  the  three  sets  of  similar  lenses  (1  and  6,  2  and  5.  and  3  and  4) 
are,  respectively,  7.825,  3.875,  and  3.820  kV.  The  envelope  parameters  of  the  output  beam 
(designated  by  subscripts  /)  at  two  locations  are  given  in  Table  3  to  elucidate  the  nature  of 
focusing,  z  =  287  mm  corresponds  to  the  end  of  the  last  lens  in  the  LEBT.  These  results 
suggest  that  the  ESQ  LEBT  system  can  deliver  a  converging,  round  beam  to  the  next  stage, 
typically  an  RFQ  in  the  linear  accelerator  section,  and  the  beam  will  be  well  matched  if  the 
RFQ  entrance  is  very  closely  (<  1  cm)  coupled  to  the  last  lens  of  the  LEBT. 

An  intriguing  feature  in  the  present  LEBT  is  that  the  combination  of  the  four  lenses 
(number  1,  2,  5,  and  6)  is  used  to  transform  the  diverging  input  beam  into  a  converging  one, 
while  the  other  two  lenses  (number  3  and  4)  act  as  a  FODO  transport  section  providing  an 
adequate  length  of  the  LEBT,  The  FODO  section  may  be  cascaded  to  obtain  any  desired 
length  of  the  LEBT. 

B.  3D  field  mapping  and  evaluation  of  fringe  fields 

The  exact  geometry  of  the  lens  system  is  used,  and  the  equipotentials  are  mapped  solving 
the  3D  Laplace  equation.  Cylindrical  electrodes  with  some  shaping  of  the  end  faces  following 
Laslett  [25]  and  Dayton  et  al.  [26]  are  used,  when  the  field  is  expected  to  be  quadrupolar 
within  about  90%  of  the  lens  aperture.  In  order  to  minimize  interference  between  two 
adjacent  lenses  [27],  a  grounded  plate  is  inserted  between  the  lenses.  Figure  2  shows  the 
equipotentials  of  the  ESQ  lens  system.  Following  the  notations  of  Matsuda  and  Wollnik 
[19],  the  potential  distribution  due  to  a  quadrupole  lens  follows  a  symmetric  behavior  as 

4>{x,y,z)  =  -<^(-x,y,2), 

4>{x,y,z)  =  -<^(i, -y,2), 

<i>{x,y,z)  =  (f>{x,y,-z), 

if  the  four  electrodes  in  a  lens  are  biased  at  alternately  positive  and  negative  potentials.  The 
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potential  distribution  can  be  expressed  as 


It  \  2\  4\i 

^{x,y,z)  =  —:^{x  -  y  ) - ~y  )  +  •••’ 


up  to  the  fourth  order,  while  the  electrostatic  field  components  are  governed  by 

E^{x,y,z)  =  -h  +  . . . , 

Ey{x,y,z)  =  k{z)y  ~^y^  - 
E^{x,y,z)  =  0-^(x2-y^)  +  ..., 

up  to  the  third  order.  The  first  term  in  the  above  expressions  represent  the  field  components 
of  a  pure  quadrupolar  field;  the  fringe  field,  which  has  a  nonlinear  behavior,  contributes  due 
to  the  second  derivative  term  in  k  for  Ex  and  Ey  and  is  governed  by  k'  for  Et.  The  function 
k{z)  is  determined  from  the  equipotentials.  This  result  in  conjunction  with  the  values  of  the 
beam  amplitude  at  various  z-locations  from  Fig.  1  is  used  to  evaluate  a  comparative  estimate 
of  the  main  field  and  the  fringe  field  which  the  beam  envelope  is  expected  to  encounter  in 
our  particular  design.  Figure  2(a),  (b),  and  (c)  show  the  field  components,  respectively,  in 
X,  y,  and  z-directions;  Fig.  3  shows  the  effective  magnitude  of  the  field  components.  It  is 
evident  that  the  contribution  of  the  correction  terms  due  to  the  fringe  field  is  very  localized 
over  a  small  region  towards  the  edges  of  each  lens  in  the  ESQ  LEBT,  and  the  main  field 
significantly  dominates.  This  suggests  that  geometrical  aberrations  may  not  be  serious  in 
the  present  ESQ  system. 

C.  Beam  dynamics  calculations  and  analysis  of  the  sources  of  emittance  growth 

The  results  from  the  above  two  steps  of  analysis  are  used  in  a  particle  simulation  code, 
PARMILA,  to  study  the  beam  dynamics  in  further  details.  This  code  essentially  pushes 
particles  using  a  chain  of  transfer  matrices,  where  the  matrix  elements  are  characterised  by 
the  factors  governing  the  motion  of  the  beam,  e.g.,  space  charge,  drift  space,  lens  elements, 
etc.  This  code  adapts  the  following  key  points. 


(a)  The  space-charge  force  is  calculated  assuming  the  beam  cis  an  ensemble  of  typically 
5000  ring-like  macroparticles.  The  macroparticles  follow  certain  distributions  in  both 
configuration  and  phase  spaces,  e.g.,  K-V,  waterbag,  semi-Gaussian,  Gaussian,  etc. 
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(b)  The  electrostatic  field  due  to  the  lens  is  decomposed  into  two  parts:  a  main  field  of 
hard-edge  type  within  the  lens  region,  and  fringe  fields  outside  this  ideal  boundary. 
In  order  to  get  a  realistic  evaluation  of  the  hard-edge  boundary  which  can  effectively 
represent  the  quadrupolar  field  component  of  a  lens,  an  equivalent  length  is  defined  as 

1 

/eff  =  —  /  K{z)dz 

Ko  Jzi 

where  kq  is  the  peak  value  of  «  in  the  flat-top  region  as  shown  in  Fig.  2,  rj  and  z-^ 
correspond  to  the  zero-crossing  locations  of  k(z)  on  the  immediate  two  opposite  ends 
of  a  lens.  A  set  of  integrals  involving  k{z),  as  prescribed  by  Matsuda  and  WoIInik  [19], 
is  used  to  model  the  force  due  to  the  fringe  field.  This  force  is  applied  at  the  ideal 
haxd-edge  boundary,  and  the  shifts  and  bends  of  the  particle  motion  due  to  fringe  fields 
are  incorporated  in  the  PARMILA  code  [22j. 

(c)  The  trajectory  of  the  macroparticles  is  followed  systematically  with  due  considerations 
of  field  distributions  over  a  cross-section  of  the  beam.  The  energy  conservation  law  is 
followed  to  tackle  this  problem 

\mvl  +  eVo  =  imt7(r)*  +  eK(r). 

where  vq  and  v{r)  denote,  respectively,  the  velocity  of  particles  on  the  axis  (this  is 
equal  to  (2eV6/nz)‘/^)  and  at  an  arbitrary  location  r,  and  Vo  and  V{r)  correspond  to 
the  electrostatic  potential  due  to  the  lens  at  the  respective  locations. 

Figure  4  shows  the  behavior  of  the  particle  distribution  at  the  output  of  the  ESQ  LEBT 
when  a  K-V  type  distribution  of  the  input  beam  particles  is  assumed.  The  distortions  in  the 
output  beam  are  due  to  nonlinear  effects,  and  this  leads  to  some  emittance  growth.  In  this 
particular  example,  the  emittance  growth  is  a  factor  of  about  1.6.  The  ellipses  composed  from 
the  Twiss  parameters  for  the  output  beam  show  a  very  similar  nature  in  the  two  orthogonal 
planes;  these  ellipses  are  drawn  using  the  effective  emittance  e  =  4e  due  to  Lapostolle  [29]. 
The  corresponding  characteristic  beam  parameters  are:  Am**  =  1.6  mm,  Fm**  =  1.3  mm, 
(An,**)'  =  -41.7  mrad,  and  (Vln**)'  =  —39.2  mrad;  a  better  matching  of  the  parameters 
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in  the  two  orthogonal  planes  may  be  fortuitous.  Various  other  types  of  distributions,  e.g., 
waterbag,  semi-Gaussian,  Gaussian,  have  been  also  considered.  The  particle  distribution  of 
the  output  beam  in  its  configuration  space  is  shown  in  Fig.  5  for  the  various  distributions.  A 
small  group  of  particles  (~  5%)  is  detached  from  the  dense  core  group  in  the  case  of  Gausian 
distributions;  these  particles  are  essentially  lost  from  the  system  as  their  location  corresponds 
to  the  boundary  of  the  electrodes.  The  rapid  blow-up  of  the  beam  for  a  non-KV  distribution 
may  be  an  artifact  of  excess  energy  in  tail  particles  of  the  distribution  function.  For  the 
sake  of  comparison  of  the  beam  behavior  in  the  two  extreme  cases,  the  output  beams  are 
shown  in  Fig.  6  for  K-V  and  Gaussian  (with  the  tail  cut  off  at  \a)  distributions  of  the  input 
beam.  The  phase-space  area  occupied  by  the  core  group  of  particles  seems  to  be  similar  in 
the  two  cases.  The  emittance  of  the  output  beam  will  practically  follow  the  dynamics  of  a 
K-V  type  input  beam  if  distribution  is  tailored  to  a  low-temperature  beam.  As  a  K-V  type 
distribution  is  commonly  used  in  the  literature  for  basic  analysis  of  beam  dynamics  and  also, 
it  seems  to  represent  well  the  beam  distribution  particularly  in  the  case  of  a  space-charge 
dominated,  high-brightness  beam  [28],  we  have  extended  our  analysis  assuming  a  K-V  type 
distribution  of  the  input  beam. 

The  evolution  of  emittance  growth,  as  the  beam  propagates  through  the  ESQ  LEBT,  is 
shown  in  Fig.  7.  It  is  noted  that  the  emittance  is  prirnarily  enhanced  at  the  second  lens  in 
the  x-plane  and  at  the  fifth  lens  in  the  y-plane.  This  result  conforms  to  the  previous  findings 
from  the  linear  beam  optics  code  when  large  beam  excursions  (Fig.  1)  have  been  noted  at 
the  aforementioned  two  lens  elements;  nonlinear  effects  act  strongly  when  the  amplitude  of 
the  beam  envelope  is  large.  In  our  simulation  studies  it  is  possible  to  diagnose  the  principal 
contributors  to  emittance  growth  and  reject  some  of  the  unwanted  group  of  particles  by 
using  appropriately  shaped  ground  plates  between  the  lenses  as  beam  scrapers,  and  thus 
the  emittance  growth  may  be  controlled.  Figure  8  shows  an  emittance  growth  by  a  factor 
of  about  1.3  only  when  a  beam  scraper  is  introduced  between  the  third  and  fourth  lenses; 
about  15%  of  the  beam  particles  is  rejected  here.  This  seems  to  be  an  easy  way  to  control 
emittance  of  the  output  beam. 
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In  an  effort  to  further  understand  the  effect  of  nonlinear  forces,  the  architecture  of  the 
simulation  code  is  made  appropriately  to  evaluate  the  contribution  of  various  factors  piece- 
wise.  i.e.,  term  by  term.  First,  the  transverse  energy  of  all  the  beam  particles  is  numerically 
maintained  at  a  constant  value  through  the  entire  LEBT.  Figure  9  shows  the  output  beam 
distribution.  It  is  evident  that  the  output  beam  is  free  from  any  major  distortions.  The 
emittance  growth  is  estimated  to  be  reduced  by  more  than  90%.  This  suggests  that  the 
chromatic  aberration  is  the  principal  contributor  to  emittance  growth  in  the  present  case. 
Turning  the  fringe-field  terms  off  in  the  simulation,  it  is  noted  that  the  emittance  growth  is 
reduced  very  insignificantly.  This  indicates  that  geometric2d  aberrations  do  not  play  a  major 
role  here. 

The  sensitivity  of  output  beam  parameters  with  variation  of  beam  voltage,  beam  cur¬ 
rent.  injection  error  has  been  studied.  Using  the  linear  beam  optics  code  it  is  noted  that 
the  beam  parameters  do  not  change  noticeably  for  a  ±  1%  change  of  quad  voltage  (from  the 
ideal  setpoint)  on  all  the  lenses  simultaneously.  Similar  insensitivity  is  noted  for  variation 
of  beam  current  (ideal  =  30  mA)  within  a  few  milliamperes.  The  effect  of  an  injection  error 
is  studied  using  PARMILA.  With  the  variation  of  the  amount  of  off-centering  of  the  input 
beam,  the  qualitative  nature  of  the  phase-space  distribution  of  the  output  beam  remains 
almost  invariant,  while  the  beam  centroid  shifts  coherently.  A  translational  off-centering  by 
1  mil  (=  0.0254  mm)  at  the  input  causes  an  off-centering  of  the  beam  by  about  3  mil  at  the 
output  of  the  LEBT;  this  also  introduces  an  angular  error  of  the  beam  centroid  by  about 
1.5  mrad.  An  on-centered  input  beam  with  an  error  in  the  injection  angle  by  1  mrad  shows 
an  off-centering  of  the  output  beam  by  about  2  mil;  the  error  in  the  angle  does  not  change 
appreciably  through  the  LEBT  channel. 

3.  Development  of  the  ESQ  LEBT  apparatus 

The  aforementioned  guidelines  are  used  to  develop  an  ESQ  LEBT  system.  Figure  10 
shows  a  schematic  of  the  LEBT  system  developed  in-house  at  Maryland.  The  principal 
elements  of  the  ESQ  LEBT  are:  (i)  six  lenses,  each  consisting  of  four  quadrupoles,  (ii) 


ground  plates  between  the  lenses,  and  (iii)  precision  ceramic  insulating  balls.  The  four 
quadrupoles  of  each  lens  are  machined  out  of  a  single  piece  of  finished  cylindrical  aluminium 
rod.  The  dimensional  accuracy  of  each  electrode  is  measured  to  be  within  ±  0.2  mil  (average) 
of  the  designed  value.  The  aluminium  ground  plates  are  2  mm  thick.  Ceramic  bails  are  used 
to  develop  an  adjustment-free,  self-aligned  lens  assembly.  Two  sets  of  precision  balls  with 
tolerance  in  sphericity  within  ±  0.025  mil  are  used.  Larger  balls  of  3/8  inch  diameter  are 
used  between  two  opposing  quadrupoles  of  adjacent  lenses;  smaller  balls  of  1/8  inch  diameter 
are  placed  between  a  quadrupole  and  its  neighboring  ground  plate.  The  ceramic  balls  are 
hidden  in  the  shadow  of  the  electrodes  in  order  that  they  remain  mostly  obscured  to  the 
field  of  the  propagating  beam.  The  entire  lens  assembly  is  mounted  in  a  cylindrical  housing, 
which  is  installed  in  a  large  vacuum  vessel. 

Voltage  hold-off  tests  have  been  conducted  on  the  lens  assembly.  With  gradual  condi¬ 
tioning  of  the  surface  of  the  electrodes  by  allowing  some  low-level  corona  discharges,  the  lens 
system  could  finally  be  driven  at  high  voltage  (up  to  about  ±  14  kV)  without  any  noticeable 
breakdown.  This  insures  practicability  of  the  present  ESQ  lens  design. 

4.  Conclusions 

This  article  addresses  the  problems  of  designing  an  efficient  low-energy  beam  transport 
system,  particularly  in  relevance  to  space-charge  dominated,  high-brightness  H~  beams.  Al¬ 
though  several  researchers  used  ESQ  LEBT  systems  in  the  past  and  showed  some  interesting 
results,  the  input  beam  parameters  were  not  as  severe  as  considered  in  the  present  case  in 
terms  of  space-charge  forces,  beam  divergence,  etc.;  further  systematic  study  of  beam  dy¬ 
namics  has  thus  been  missing.  An  emphasis  is  given  here  on  simulation  studies  to  examine 
some  beam  physics  issues  critically.  This  approach  lays  a  foundation  to  develop  an  efficient 
LEBT  system;  especially,  the  role  of  the  various  nonlinear  forces  on  emittance  dilution  is 
studied  in  depth. 

The  lens  assembly  is  designed  following  several  iterations  of  the  numerical  scheme  so 
that  the  lens  geometry  is  optimized  as  much  as  possible.  Typical  beam  characteristics  of  a 
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high-brightness  H“  source,  Penning- Dudnikov  type  (as  used  in  the  BEAR  experiment  (4]), 
are  used  in  the  analysis  with  an  aim  of  transporting  a  space-charge  dominated,  diverging 
beam  from  the  source  over  a  length  of  about  30  cm.  and  yield  a  converging  beam  at  the  end 
without  significantly  deteriorating  the  intrinsic  emittance  of  the  beam.  The  beam  dynamics 
is  studied  in  detail  using  a  particle  simulation  code,  which  enables  us  to  study  the  influence  of 
various  factors  in  determining  the  emittance  growth,  e.g.,  nonlinear  effects  due  to  chromatic 
aberrations  and  fringe  fields,  various  types  of  distributions  of  the  beam,  sensitivity  of  various 
control  parameters,  etc. 

The  mzdn  contributor  to  emittance  growth  in  the  ESQ  LEBT  is  found  to  be  chromatic 
aberrations.  It  is  shown  in  the  simulations  that  particles  responsible  for  enhancing  the 
emittance  growth  can  be  identified  and  subsequently,  they  are  rejected  by  carefully  using  two 
intermediate  ground  plates  as  beam  scrapers.  Thus  the  emittance  growth  can  be  controlled 
without  sacrificing  the  beam  current  appreciably.  In  the  case  of  the  H"  beam  from  the 
Penning-Dudnikov  type  source  the  nominal  emittance  growth  through  the  ESQ  LEBT  system 
is  estimated  to  be  a  factor  of  about  1.6,  while  this  factor  is  reduced  to  about  1.3  by  rejecting 
only  15%  of  the  beam  particles.  Such  a  modest  value  of  emittance  growth  attributes  merit 
to  the  particular  design  of  the  ESQ  LEBT. 

The  ESQ  LEBT  system  developed  in-house  at  Maryland  own  some  uniqueness  in  its 
mechanical  assembly.  The  entire  system  is  self-aligned  and  it  is  very  rugged  mechauiically. 
Each  component  is  fabricated  with  a  dimensional  tolerance  close  to  the  designed  value. 
Results  on  voltage  hold-off  tests  of  such  a  compact  system  are  encouraging;  no  noticeable 
breakdown  events  occur  up  to  a  voltage  level  of  about  ±  14  kV. 

Beam  transport  experiments  are  planned  using  the  ESQ  LEBT.  Details  of  the  experi¬ 
mental  work  and  further  studies  on  beam  dynamics  will  be  reported  elsewhere. 
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Table  1.  Estimated  H“  beam  parameters 


Beam  current  4  30  mA 

Beam  voltage  Vj  35  kV 


u/c  =  l3 

Generalized  beam  perveance  K{=  2Ii,l 
Initial  rms  normalized  emittance 
Initial  beam  radius  a 

Initial  divergence  of  beam  envelope  at  r  =  a 


8.6  X  10-^ 

3  X  10-^ 

0.069  IT  mm-mrad 
1  mm 
20  mrad 


Table  2.  ESQ  lens  parameters. 

Lens  Number 


1  &6 

2  &  5 

3  &  4 

(mm) 

8.00 

12.00 

12.00 

1  (mm) 

15.00 

59.00 

47.00 

L  (mm) 

6.00 

6.00 

6.00 

Table  3.  Beam  parameters  at  the  output  end  of  the  LEBT. 

z  (mm)  Xf  (mm)  Yj  (mm)  Xj  (mrad)  Yj  (mrad) 

293.0  1.39  1.37  -53.1  -53.7 

302.0  1.03  1.00  -24.0  -24.3 


Figure  Captions: 

Fig.  1.  X  and  Y  are  the  amplitudes  of  the  beam  envelope  in  .x  and  y  directions;  z  is  the 
direction  of  propagation.  /Cy  =  — #Cr.  Note  that  (and  Y)  =  12  mm  corresponds  to 
the  mciximum  aperture  of  the  electrodes. 

Fig.  2.  3D  LAPLACE  results  of  the  field  components:  (a)  Ex  ,  (b)  Ey  ,  and  (c)  Ex  .  Solid 
line:  pure  quadrupole  field;  dashed  line:  fringe  field. 

Fig.  3.  Effective  magnitude  of  the  fields.  Solid  line:  quadrupole  field;  dfished  line:  fringe  field. 

Fig.  4.  Particle  distribution  in  the  output  beam  (right  column)  of  the  ESQ  LEBT  for  a  K-V 
type  input  beam  (left  column).  Note  that  X  (and  Y)  =  8  nun  corresponds  to  the 
aperture  of  the  first  and  last  lenses. 

Fig.  5.  Configuration-space  geometry  of  the  output  beam  for  various  input  distributions:  (a) 
K-V,  (b)  waterbag,  (c)  semi- Gaussian,  and  (d)  Gaussian. 

Fig.  6.  Comparison  of  phase-space  distribution  of  the  particles  in  the  output  beam  for  K-V 
(left  column)  and  Gaussian  (right  column)  input  distributions. 

Fig.  7.  Spatial  evolution  of  rms  emittance  through  the  ESQ  LEBT.  The  emittance  values  are 
plotted  at  the  end  of  each  lens  in  the  LEBT. 

Fig.  8.  Phase-space  distribution  of  particles  in  the  output  beam  when  15%  of  the  beam  is 
intercepted  by  two  beam  scrapers. 

Fig.  9.  Particle  distribution  in  the  output  beam  when  all  the  particles  are  assumed  to  be 
transported  with  the  same  energy  (35  kV). 

Fig.  10.  Schematic  layout  of  the  ESQ  LEBT. 
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Abstract 

A  moment  method  technique  for  solving 
Laplace's  equation  is  presented.  The  technique  is  then 
extended  to  Poisson's  equation  to  include  space  charge. 
The  procedure  is  implemented  on  a  PC  and  applied  to 
the  cases  of  an  electrostatic  quadrupole  (ESQ)  lens  and 
an  ellipsoidal  bunch  in  a  grounded  pipe. 

I.  INTRODUCTION 

Currently  we  are  engaged  in  the  design  and 
development  of  a  Low  Energy  Beam  Transport  (LEBT) 
section  for  H'  ion  beams  [I].  We  have  chosen  to 
employ  electrostatic  lenses  in  this  design,  since  for  low 
ion  velocities  magnetic  lenses  fail  to  provide  enough 
focusing  for  intense  beams  while  gas  focusing  is 
intrinsically  stochastic.  The  current  prototype  consists 
of  6  ESQ  lenses.  In  order  to  model  the  action  of  such 
a  lens  it  is  necessary  to  solve  Laplace's  equation  for  the 
particular  lens  geometry.  Once  the  lens  is  characterized 
electrically,  the  information  may  be  used  in  other 
simulation  tools  to  aid  in  design. 

In  general,  numerical  methods  must  be  utilized 
to  solve  Laplace's  equation.  We  present  a  technique 
which  is  fully  three  dimensional  yet  is  efficient  enough 
for  implementation  on  a  PC.  The  efficiency  of  the 
technique  arises  from  the  fact  that  it  is  based  on  an 
integral  formulation  rather  than  the  more  common 
differential  form.  Instead  of  solving  for  electrostatic 
potential  directly,  we  solve  for  the  surface  charge 
density  on  conducting  bodies.  This  results  in  a  reduced 
dimensionality  of  the  problem  domain.  The  integral 
formulation  also  readily  extends  itself  to  Poisson’s 
equation.  Thus,  we  can  model  lenses  in  the  presence  of 
charge  distributions.  Also,  since  we  know  the  surface 
charges,  we  may  evaluate  capacitances  between  various 
lens  elements. 

II.  NUMERICAL  TECHNIQUE 

A.  Laplace's  Equation 

Letting  (|)  denote  electrostatic  potential,  the 
problem  is  usually  seen  in  the  mathematical  form 


=  0  V;teQ, 

<Mx)  Vxer. 

Here  Q  is  the  3D  region  of  interest  and  1'  is  its  boundary  (i.e. 
r=5Q).  The  function  f  represents  the  given  boundary  values 
and  constitutes  the  data  of  the  problem.  Usually  Q  represents 
the  beam  line  and  F  is  the  surface  of  a  focusing  lens,  thus  f 
would  be  the  lens  voltage.  A  finite  differencing  method  would 
typically  attack  this  problem  directly.  However,  we  prefer  to 
work  with  an  integral  representation  of  the  problem  rather  than 
the  differential  form  [2|. 

/X)  =  /C(x,Uo(«d-5,  (2) 

r 

where 


G(x,0  = 


1 

4it  |x  -  ^  1 


(3) 


is  the  free  spaee  Green's  function  for  Poisson's  equation.  The 
function  a  is  introduced  as  the  new  unknown  for  the  problem. 
It  is  recognized  as  the  surface  charge  density  on  the  boundary 
r.  Once  o  is  known,  4>  may  be  recovered  \  ia 


4»(x)  =  fG(x,^)a(OdH  VxeQ.  (4) 

r 

Note  that  the  dimensionality  of  the  problem  has  been  reduced. 
In  (1),  must  be  solved  on  fj,  a  3D  subset  of  E*,  while  (2)  is 
defined  only  on  the  2D  manifold  F. 

We  employ  the  method  of  moments  to  solve  (2),  the 
details  of  which  are  presented  in  |3].  Qualitatively,  the 
technique  is  vary  similar  to  the  representation  of  a  quantum 
mechanical  operator  in  matrix  form.  We  chimsc  a  set  of 
expansion  functions  {u,}-  on  F  which  is  used  to  approximate 
CT.  that  is  o(x)=Za,u,ix),  for  some  set  of  a.e  R.  Another  set  of 
functions  called  weighting  (or  testing)  functions  {v„;  is  also 
selected.  After  expanding  (2)  in  |u,J  we  take  the  inner 
product  with  each  of  the  v^'s.  We  end  up  with  a  scries  of 
linear  equations  where  the  a,'s  are  the  unknowns.  This  system 
may  be  solved  by  standard  matrix  methods.  We  apply  a 
conjugate  gradient  algorithm  to  this  end  |4|.  this  is  an  iterative 
method  which  seems  to  provide  fastest  convergence. 

For  our  moment  method  we  chose  for  {ii.J  a  set  of 
piecewise  constant  functions,  constant  mcr  the  face  of  a 
triangle.  Specifically,  F  is  triangulated  (appniximated  by 
triangles,  for  example  sec  figure  I )  and  a  is  as.sumcd  constant 
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over  each  triangle.  This  selection  results  in  a  tinite 
element  representation  to  (2).  For  the  weighting 
functions  wc  selected  Dirac  delta  functions  located  at  the 
centroid  of  each  triangle  (this  is  known  as  point- 
matching).  This  allows  fastest  evaluation  of  the  inner 
products  and  yields  good  results  as  long  as  the  triangles 
are  sufficiently  regular. 

B.  Extension  to  Poisson's  Equation 

If  we  wish  to  model  a  charge  distribution  p  in 
the  presence  of  our  boundaiy  T,  it  is  convenient  to 
exploit  the  linearity  of  the  integral  operator  in  (2).  That 
is  the  potential  at  the  boundary  must  be  the  sum  of  that 
due  to  both  a  and  p. 

Ax)  =  /c?(x.5)o(5)d^5  *  (5) 

r  Q  ^0 

or 

Ax)  -  <t),(jc)  =  (6) 

r 

where  (tip  is  the  free  space  potential  due  to  the  charge 
distribution  p,  given  by  the  second  integral  in  (5). 
Equation  (6)  is  similar  in  form  to  (2)  and  may  be  solved 
by  applying  the  method  of  moments  as  before  to  the 
boundary  data  f-<j)p  .  Therefore,  it  is  only  necessary  to 
determine  itip,  the  free  space  potential  due  to  p,  in  order 
to  apply  the  method  to  Poisson's  equation.  This  may  be 
done  numerically  or  analytically  (if  available). 

III.  APPLICATIONS 


The  preceding  technique  was  implemented 
using  Borland  C++  3.1,  The  platform  was  an  t486  PC 
operating  at  33  MHz  and  running  Windows  3.1 
operating  system.  All  examples  where  run  in  double 
precision  arithmetic. 


-4.  ESQ  Lens 

A  3D  potential  problem  is  the  modeling  of  an 
electrostatic  quadrupole  lens.  Figure  I  shows  the  computer 
model  of  an  ESQ  lens  similar  to  the  type  used  in  [1].  It  is 
formed  from  4  cigar-shapied  electrodes,  the  beam  would  enter 
from  the  left.  Each  electrode  is  59  mm  long  and  has  a  radius 
of  12  mm.  The  aperture  of  the  entire  lens  is  10.5  mm.  Two 
grounding  shunts  are  located  at  z=+3 1  mm  (they  are  not 


Figure  2:  ESQ  Lens  Focusing  Function 

shown  in  figure  I  to  avoid  clutter)  which  provide  isolation 
from  adjacent  lenses. 

The  single  particle  focusing  effect  (  the  kappa 
function  ic(z) )  from  such  a  lenses  can  be  determined  from  the 
derivatives  dEydx  and  dE/dy  on  axis.  Figure  2  shows  the 
computed  data  for  the  y-plane  for  the  case  in  which  the  x- 
plane  electrodes  are  driven  to  IV  and  the  y-planc  electrodes 
are  held  at  TV.  The  grounding  shunt  at  cither  end  of  the  lens 
cause  the  rapid  decay  in  dE/dy. 

B.  Ellipsoidal  Bunch  in  a  Pipe 


Figure  J  2. 51  Ellipsoid  in  Pipe 


We  can  use  the  Poisson  extension  to  simulate 
a  uniform  charge  density  ellipsoid  in  a  conducting 
cylinder.  This  situation  is  useful  in  modeling  cold 
bunched  beams  propagating  through  a  beam  pipe. 
There  exists  an  analytic  solution  for  the  potential  of  such 
an  ellipsoid  in  free  space  [^.  Thus,  it  is  only  necessary 
to  model  the  pipe  (surface  charge)  numerically. 


Figure  4;  10:1  Ellipsoid  in  Pipe 


Figure  3  shows  the  axial  potential  and  z 
component  of  electric  field  for  the  case  of  an  axis- 
symmetric  ellipsoid  with  major  axis  2.3  cm  and  minor 
axis  1  cm  inside  a  pipe  of  radius  2  cm.  In  this  case  the 
fields  are  still  relatively  linear,  the  image  effects  from 
the  pipe  arc  slight.  However  in  figure  4  we  see  a  very 
nonlinear  field  for  the  case  of  an  ellipsoid  with  major 
axis  10  cm  and  minor  axis  1  cm.  We  find  that  after  the 
bunch  length  becomes  comparable  to  the  pipe  radius  the 
image  effects  play  an  increasing  role.  The  total  charge 
in  both  cases  is  10'"  C. 

The  above  simulations  can  be  used  to  determine 
the  so  called  "g-factor"  for  bunched  beams  in  cylindrical 
pipes.  A  detailed  discussion  of  these  results  can  be 
found  in  [§]. 

IV.  CONCLUSION 


When  dealing  with  the  situation  of  conductors  in  a 
vacuum  the  moment  technique  is  in  general  quite  successful. 
However,  if  many  dielectrics  arc  present  it  is  probably  best  to 
use  finite  differencing.  Also,  when  it  is  necessary  to  know 
the  fields  over  a  large  set  of  points,  say  when  doing  many 
particle  simulations,  it  is  probably  best  to  use  finite 
differencing.  Computing  the  potential  is  a  moderately 
expensive  process,  since  we  must  evaluate  (4)  at  each  point, 
while  finite  differencing  solves  for  the  potential  directly. 

For  the  situations  discussed  the  technique  is  well 
suited.  We  need  full  3D  solutions,  yet  only  for  the  case  of 
conductors  in  a  vacuum.  Also,  we  are  only  concerned  with 
the  solution  data  along  the  bcamline  axis.  Therefore,  the 
number  of  data  points  to  compute  is  a  minimum. 
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The  method  of  moment  technique  has  several 
advantages  and  disadvantages.  The  overall  advantage  of 
the  technique  stems  from  the  fact  that  only  surfaces  are 
considered,  rather  than  3D  regions.  Hence,  it  is  a  good 
method  to  model  complicated  or  otherwise  arbitrary 
geometries.  For  the  same  system  order,  we  get  a  higher 
boundary  resolution  as  compared  to  finite  differencing 
on  a  grid.  Also,  since  the  surface  charges  arc  solved  for 
it  is  possible  to  calculate  the  fields  anywhere  in  space, 
without  interpolation.  This  fact  allows  us  to  apply  the 
method  to  unbounded  situations. 


